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FOREWORD 


Thl* Is a combined report for the two project* at the Maa*achu*ette Insti¬ 
tute of Technology which are spon.ored by the Office of Naval Re.earch under Contract 

N5orl60. 


Project on Machine Method* of Computation and Numerical Analy*l* 

This Project Is an outgrowth of the activities of the Institute Committee 
on Machine Methods of Computation, established In November 1950. The purpose of the 
Project Is (1) to Integrate the efforts of all the departments and groups at M.I.T. 
who are working with modern computing machines and theli applications, and (2) to 
train men in the use of these machines for computation and numerical analysis. 

People from several departments of the Institute are taking part In the 
project. In the Appendix will be found a list of the personnel active In this pro¬ 
gram. 

Project Whirlwind 

This Project makes use of the facilities of the Digital Computer Laboratory. 
The principal objective of the Project Is the application of an electronic digital 
computer of large capacity and very high speed (Whirlwind I) to problems in mathe¬ 
matics, science, engineering, simulation, and control. 


The Whirlwind I Computer 

Whirlwind I Is of the high-speed electronic digital type, In which quanti¬ 
ties are represented as discrete numbers, and complex problems are solved by the re¬ 
peated use of fundamental arithmetic and logical (i.e., control or selection) opera¬ 
tions. Computations are executed by fractional-microsecond pulses In electronic 
circuits, of which the principal ones are (1) the flip-flop, a circuit containing two 
vacuum tubes so connected that one tube or the other Is conducting, but not both; (2) 
the gate or coincidence circuit; (3) the magnetic-core memory, In which binary digits 
are stored as one of two directions of magnetic flux within ferro-magnetlc cores. 

Whirlwind I uses numbers of 16 binary digits (equivalent to about 5 deci 
mal digits). This length was selected to limit the machine to a practical size, but 
It permits the computation of many simulation problems. Calculations requiring 
greater number length are handled by the use of multiple-length numbers. Rapid- 
access magnetic-core memory has a capacity of 32,768 binary digits. Present speed 
of the computer Is 40,000 single-address operations per second, equivalent to about 
20,000 multiplications per second. 
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PART I 


COMMITTE ON MACHINE METHODS OF COMI 


UTATI0N AND NUMERICAL ANALYSIS 


General Comments 


machine coStl^Um^T^ '7^ »« *h. 

numbers of the staff of the present protect will 7 n " r ‘ n ’’ the spring many 

TnLTu LT er ' r i p rlWind WU1 ^ WlthdraWn -e r;”, 0 e f .rch e .n e d nter 

will be turned over to Project Lincoln for Its Tull time use. 

Progress toLejT.mM"* Project will terminate and these 

Center will tike un\h T v n6W 861 1 rogre8s Sports. from the Computation 

nter will take up the task of reporting on research and training In the use of 

igital computers at MIT. All those receiving this report, who wish to have copies 

who is "““be Pr ° BreSS lte P°rts, should write to the undersign®,I 

wno is to be Director of the Computation Center 

,.„ t A “ Pa f t ° f th * tr-,nBfer ot activities, no part time ONK research asaia 

lns!L7 re r aPP ^ thi ' aC “ demlc >' oar ' IBM research assistants were appointed 
notes n t ^ 0n8e< ’ Uent ^' the fl r»l half of this Report, which usually contained 
notes on the research of the 0N11 Assistants, contains final reports from last year's 
men, who finished last summer. The next report of thi. series and al“ Z\l 
naming two will only contain what is the second half of this report, a report on 
the research activities of Whirlwind. The activities of the IBM research a.sla- 
tants and of the staff of the new Computation Center, will be reported In the first 
issue of the new series of reports, which will be Issued early in 1957. 


Philip M. Morse 

Director, MIT Computation Center 
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2.1 Ind«x to Report! 

Terminating Report: Phyelcm Department 
Terminating Report: Mathematic! Department 
Aerodynamic Equation* 

Spherical Wave Propagation. 

Thin Airfoil Tlieory: Perturbation Solution*. 

Non-Linear Boundary Value Probleaa 
Spherical Wave Propagation. 

The Problea of Buckling of a Hyperbolic Paraboloidal Shell 
Loaded by It* Own Weight. 

Spherical Wave Propagation 

Aerodynamic Theory. 

Application of Information Theory to Machine Computation*, 
quadratic Theory of Bending for a Circular Plate. 

Machine Solution of the Diffusion Equation. 

Calculation of Numbers of Structures of Relation* on Finite Bets. 

The Number Distribution of Electrons and Photons in Cascade 

The Basic Problem of Numerical Analysis Expresssed in the Language 
of Computing Machines 

Multiple Prediction 

Deterministic Processes and Stochastic Processes. 
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- 2 Terminating Reports 
PHYSICS DEPARTMENT 


contract N5ori60 reveals a remarkablv well .a . 1 pnysics under 

forts were primarily concerned"ithor^Lr? pattern ° f * ro *‘ h - 

cess which required a period of about one year" "“'h 1 "! ‘’“■Put*! Ion , a pro 

was sufficient understanding of machine processes so ^hat"fa? 1 th ‘ S ‘ h * re 

could be attempted, and, indeed, many were carried thr * f “ lrly com > ,lex problems 
Not a few of the early problems led to the rr ough to successful conclusions 

the end of the project made cornu,!,"^tlon of subroutines which, toward 
que!tion during the earlier days. 003 poSSible " hich would have been out of the 




wnu ui 


and it rapidly became the center of Interest ltoler"'‘ Ch , lne n 4< ‘ a ' ,sUable 

machine computation; mathematical analysis and h h ’ * aU rt8earch led to 

ficient for the completion of c«t“S pwjecE. C °" putat ‘ on —*!■.. suf 

group was JUly 195l> att8ntl °" of the physic 

evaluation o( . “ ,1 f , pr ° ble " a tt8 matr ‘ X *»*««<». the 

equation to obtain llT 3 iTl'Zi'£,1^**'?™'* ° f a 

Quarterly Progress ReLrtsT “Th " * ra ? dUflractlon *° rk (8 « 8 

. . • * 3 ), Thi! work was all begun on the I R M 

machines hut was transferred to Whirlwind I when that machine became svall.hle late 


k— ~vf 

‘rrp r*jr.!2 rrsr 

* a" f: 13) and J - L ‘ Uret8ky5 CQPR 5 9) carried out variational calculations 

case He uV Vh‘ tterl "f of slow electrons by atomic Hydrogen and. in the latter 
eml. L ,US " Te,8kin <QPR 16) p8rfor, * ed 8 “ch a calculat Ion for 

(OPR 7 81 [IT yK ° n ! ' so “ ewhat supplementary Investigation by L\ Mower 7 

nC.I'.hen^!o!rr, , e*t.bll.hlng the accur.cy of cert.in van.,Iona. t.ch- 

d v d , ° elytron atom scattering problems. Mower s calculation, -ere 

by hand, all others were carried out on Whirlwind I. 

w*. B°“ nd 8tBt8 Problems -ere not neglected. One of the earliest calculations 

Whlrlwlnd°to i t"""! 1 ' K <0P * 5- 8> * h0 flt a Mor8e Pot®«>tlsl to the deuteron using 
H vn 9 ? ( nte Erate the Schroedlnger equation and Iterate for the eigenvalue 

to „,“h .1 . J 3> lBprovod the t«chnlques of Morse. Young and Hurwlt.lO in order 

ltd. ! “^ulations of single particle eigenfunctions for many electron 

, A - Tubia (QPR 14 16) mechanised the procedure so that the relevant quan 
u«dl!°! ! be * enerated by Whirlwind on a demand basis M Rotenberg 12 (QPR 15 18) 
used Whirlwind to carry out a Ilartree Fock calculation for a heavy nucleus 
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Two investigations in solid state physics were accomplished under the 
contract. F. J. Corbato 13 (QPR 12 19) undertook the machine calculation of energy 
bands in graphite Professor PM. Morse (OPR 16) developed a theory for obtaining 
electronic wave functions in a crystal lattice. Calculations relative to the latter 
investigation were initiated by H Fields and Z. Pried (OPR 18). Corbato s work, in 
particular, has led to the development of many subroutines which have been useful to 
other programmers (for examples, see the QPR s referred to above). 

One physics assistant chose to do his work in the field of operations re¬ 
search. J. D. C. Little 14 (QPR 7, 8) studied the problem of optimizing the water 
use in a large hydro-electric system. Again, Whirlwind was used to carry out the 
computations . 

Two long term table making projects were initiated by the physics group. 
Sartori, Tomkin, Tubis and t/retsky were to develop a Whirlwind Program for producing 
a tabic of coulomb wave functions (QPR 2-19). Corbato and Little wero to do the 
same for spheroidal wave functions (QPR 3 11). The first project is still In pro¬ 
gress, and the second one has recently achieved publication 1 ^. 

The development of a program to generate the expansion coefficients and 
eigenvalues 1 ^ for the spheroidal wave functions made it possible to carry out a 
subsequent research problem in quantum physics. J. L. Uretsky used the basic 
routine of Corbato and Little in order to develop a technique for computing the 
bound states and zero energy scattering cross sections for a finite square well 
potential of spheroidal symmetry 16 (QPR 12-17). 

The survey Just given is not all inclusive. As in any research program 
there were many avenues of Investigation which were abandoned for one reason or 
another. Some of these were: 

Programming of the Knight a Tour ' on Whirlwind (QPR 2) 

Scattering of Electrons from a Diatomic Molecule (QPR 5, 6) 

Neutron Proton Interactions (QPR 5, 6) 

Whirlwind Study of Group Dynamics (QPR 5) 

Photoproduction of Mesons from Deuterons (QPR 8, 11) 

It seems appropriate to remark that many members of the faculty of the 
M.I.T. Physics Department played indispensable roles in initiating and guiding the 
research carried out under this contract. In particular, Professors P. M. Morse 
and H. Feshbach are, in one way or another, associated with nearly all of the in¬ 
vestigations discussed above. 
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5 iurXSL SJSS? a= **** 

6. A. Teak In, Ph.D. Thesis Physic. (title not available). 

7 ' Pot£;::; ***** * sc.tt.nng fr«. s ,. uc 

t ntl»l- M.J.T. Physics Department. Ph.D. Thesis (19537 - ‘ - - — 

Physlcni Review 89, 947 (1953); 99. 1065 (1955). 

8 w'i C £" t fh liC ,’ - ^l PO,1Polo > lc *L Triplet Neutron-Proton Potential. 

M.I.T. Physics Department, Master s Thesis ( 1953 ). -- 

9 ' N I Y T 1 "^^7^ ? n? in - C ^ £5i — ^ ,n - ,ltlon Probabilities of Light Atoms 
M.I.T. Physics Department, Ph.D. Thesia (1954)--~- 

Physical Review 100, 1148 (1955). 

10. P. M. Morse, L. A. Young, and E. S. Haurwitx, Physical Review 48 . 948 (1935). 

11. A. Tubis, Physical Review 102 , 1049 (1956). 

12 I'.TTT A Self-Consistent Determination of the Nuclear Radius. 

S' T' C ^ ttee on r Mr^' Z 0 ™ e ‘ iS PhF.i^iT iSvIWTob. 439 (1955) 

M.I.T. Committee on Machine Computation Technical Report No. 6, S^fember 1955 

13. F. J. Corbett?, J. Chem. Physics 24, 452 (1956). 

14 ' — — orilg - e - W,ter iS » "ydro-elec tr!c System, 

)^nca Ho’ 2 mw" 1 ' Ph ' D ' 11,0,18 fl955) - ^ Soc - 

15 ' V’ “° r8Cl L ‘ J ' ChU ’ J ' “• C ' Wttle. and F. J. Corbato. 

Spheroidal Wave Functions . Wiley (1956). 

16- i\ L I U " t8ky ’ 525S Properties of the Spheroidal Square-Well Potential 
M.I.T. Physics Department,Ph D, Thesis (lgs'eTT - Ph7s"lcil Rev 1^7 to bi 
published. 


Jack L. Uretsky 


MATHEMATICS DEPARTMENT 
Introduction 


Foo tnotes 

Quarterly Progress Reports of the M.I.T. Coaimlttee on Machine Computation and 
Numerical Analysis will be referred to by the initials 'QPR' 1 followed by the 
relevant report numbers. Report number 1 was dated September 15, 1951. 

H. Feshbach and S, I. Kublnow, Physical Review 88, 481 (1952); 96, 218 (1954)- 
98, 183 (1953). — ' 

L. Sartort, Ph. D. Thesis - Physics - (title not available). 

M. C. Newsteln, A Va r iatio nal Treatment of Atomic Scattering. 

M.I.T. Physics Department. Ph.D. Thesis <1955). M.I.T. cSSTttee on Machine 
Computation, Technical Report No. 4, February 1955. 


hv »h 18 8 ternlnatin 8 report for work done under DIC project number 6915 

Dy tne mathematics members during the period July 1, 1951 to June 15, 1956 The 
research assistants and associates selected from the Mathematics Department and 
constituting all the mathematics members except faculty - that have worked under 
ponsorship of the Committee at any time during the above period have been the 
following: 

Aronson, Donald G. 

Benlngton, Herbert 
Burrows, James W. 

Caldwell, William V. 

Fox, Phyllis A, 

Glantz, Herbert 
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Handy, Benjamin P. 

Hershenov, Joseph 
Hicks, Noel J. 

Hsu, Pao Tau 
Klein, Paul 
Lister, Mary 
Mcllroy, M. Douglas 
Phipps, Philip L. 

Ralston, Anthony 
Rankin, Bayard 
Schlesinger, James 
Sogel, Lee 

Welnitschke, Hubertus J. 

The main body of this report is made up of abstracts of selected projects that were 
inaugurated and carried out by ten of the above members. "Hie selection of projects 
to be abstracted was made largely on the basis of how completely the work was re¬ 
ported and how well the project was completed or led the researchers to clearly 
stated results This selection has left out any mention of a fair amount of work 
that in the course of the projects warranted progress reports but did not seem to 
lead to any conclusive results. Some of the work done under the committee is too 
now to report in conclusive form and has not been mentioned for that reason. The 
projects that have been abstracted encompass the subjects of aerodynamic theory, 
fluid dynamics, elasticity, differential equations, group theory, cascade shower 
theory, stochastic processes, and the logic and programming of digital computers. 

The abstracts themselves have tried to outline the particular problem at 
hand, the mathematical structure of the problem, the techniques used to solve the 
problem, and the form of the results Mention has been made of the type of digital 
computer used, if any. After the title of each abstract some hyphenated numbers 
identify the references within the progress reports where the work was reported. 

For example 1-11 refers to progress report no. I, page 11. Further references in 
the form of final reports, technical reports or published papers are numbered in 
square brackets at the end of each abstract and accumulated at the end of this 
terminating report. 
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,lo»ai q varl a blo a p variables'*"^ 7*2**^ 

pairs of equation, that were ultinately n .olved eSral ° qUatl0nS and 11 was the9e 


Fox * gP her * c »l Hive Pro p agation . 2-12, 3-6, 4-9. 


equations. A the n inteKraUon C of e eoLtl^nI 1Cal ln * egratlon of hyperbolic differential 
adiabatic expansioni st TIZZZ »» 

expressing the equations in fm^ carried out on Whirlwind I by 

'“ ,r «~- 

... **• °<.................. 

1 - 9 > . - C* 

-g- + d~ a _ ^ dt 

Along dr = (u - c)dt they take the form 

iLf. - dr“-_ 3 - J-U. dt 

/o C V 

15 ‘ hBSe latt ? r equatlona ^Pressed in difference form that were solved Im- 
Wh7e CO “ PUtatlonal difficulties arose in regions where terms behaved badly ur 

where P ° nV ® rB ° nce * as P°° r - Such regions were at the center of the sphere and 
where shock waves developed. K 


Abs tracts 

P. T. Hsu, A erodynamic Equation s. 1-11, 29, 3-7, 4-9. 


This work was later carried on by 
p. Fox it led to part of the first technical 
Ph.D. Thesis in mathematics. 


A. Kalston. Under the authorship of 
report of the project [ lj and to a 


The general problem of this project was to obtain approximate solutions ' Fox ’ Thlp Air foil Theory : P erturbation Solutions . 6-16, 7-17, 8-20. 

of the two dimensional lineai Integral equations of lifting surface theory for 

steady Incompressible flow. The problem was confined to the study of a rectangular A P er turbation method was used to solve some problems of fluid flow past 

surface and had interpretation for 2 full Bpan aileron. Specifically, the pairs “ th * n alrfoil • The Several method that effects a linearisation of the problem 

of quantities desired were the lift and moment, or the lift and hinge-moment. Prof. *“* ia gm> ved with a technique described by M. J. Llghthlll. Solutions were obtained 

E. ftelssner suggested the method of solution. The numerical technique required was mainly for the simple two-dimensional, irrational, incompressible flow but some 

largely numerical integration. Calculations were carried out at first by desk cal- results for the more interesting compressible case were obtained. The’ehape of the 

culator and later by digital computer. Numerical results were obtained and compared airfoil was chosen as an ellipse though some consideration was given to the possible 

with those of Weissinger, another worker In the field. application of the Lighthlll technique to hyperbolic types. 


Mathematically speaking, the basic equation was the following which was 
to be solved for p when w and k were given: 



The equations Involved are, for the velocity components u and v. 


u 

X 


+ v 


y 


o 


U -V * 0 

y * 
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The Lighthlll technique involves a coordinate transformation based on expanding the 
independent as well as the dependent variables in terms of a small parameter. The 
consecutive functional coefficients of the powers of the parameter are adjusted so 
as to suppress singularities arising in the solution of the problem. This Involves 
a transformation of the coordinate system. The chief difficulties were encountered 
in the convergence of the modified perturbation series near the leading and trailing 
edges of the airfoil. 

This work led to related techniques In the study of plane wave propagation 
|2j and to part of the technical report [ lj and thesis already mentioned. 

N. J. Hicks, Non-Linear Boundary Value Problems . 4-8 # 6-9, 7-14, 8-16. 

Some series solutions for boundary value problems of ordinary non-linear 
differential equations were carried out on Whirlwind 1. The equations of interest 
appear in the theory of finite deformations of thin elastic circular plates. These 
equations were solved by expressing the unknowns in Taylor series, obtaining recur¬ 
rence relations for the coefficients in these series, and then solving these recur¬ 
rence relations. In order for this method of solution to go through, it was neces¬ 
sary for some boundary conditions to be replaced by additional initial values 
obtainel through non-linear interpolation. The use of an I.B.M. card machine or a 
differential analyser for this problem was also considered. 

The specific equations of interest to be solved for f(x) and g(x) were 

f (x) l-f'd) - KjU) = u 

x Jr x 


g (x) + ; g'(x) - ^ig(x) = 5 * - (x) 


where u is constant. The fora of Taylor series used for substitution in these e- 
quations were 

* ‘ 5 b n <- 1 > n 


f(x) 


nO “n (X - X) 


A. Ralston, Spherical Wave Propagatio n. 6-11, 7-17, 8-18, 10-10, 11-14, 12-37. 

The work discussed above under the same title and initiated by P. Pox 
was completed. That Is, the numerical solution of the two partial differential 
equations governing the propagation of spherically symmetrical waves was carried 
to the point whore the state of the system was very close to atmospheric equili¬ 
brium and was quickly approaching it. The initial density distribution for which 
complete results were obtained was,* = 1 + 2e-«P 2 . The results showed that the 
shock, which it was thought might form, does in fact not form. The motion was 
similar qualitatively to that found by J. J. Unwin in investigating the same pro¬ 
blem. However, the compression at the center due to the recoil wave was a factor 
of three times less than Unwin's. 

In order to obtain complete results with the initial density distribution 
specified above, it was necessary to introduce near the t-axis an interpolation - 
extrapolation procedure to modify the method of two-characteristics previously used 
by Cox. This modification was then incorporated in a comprehensive program As a 

[nterLi t of a »h P ° B ? ible ^ ^ Pr °* ra “ to c °“P ute »s »«ny points as desired in the 
interior of the plane and on the t-axis for any r-axis spacing. Two integral, re¬ 
lating to mass and density were evaluated numerically in order to provide a physi¬ 
cally meaningful check on the accuracy of the results. A mathematical check on the 
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sitlvlty Of the solution to the initial r-axls spacing was also made. 

for any initial density district ^“ p " hBnslve Program shoulg have given results 
latter distribution J. not forthcominv h ” 1" 1+5 ®~ ' Kesult, for this 

A careful investlgation intr »» owever » because errors became too large, 

lions were stopped It appeared that '^ ° theS ® orrors * as made. Before calcula- 

not formed with the prevTu.don.Uy '' * h ° U ‘ b “ **- 

project faf 1 * ’° rk ” S r ° POrt0d I" the second technical report of the 

L K We^: ^ ‘^^^^Ung.^^perbol tc Paraboloidal Shell Loa ded. 

load which "causes **t he ~ ^ 8 " aU * 8t <“«tributed 

moment free suppor „d To.l ^ hvp8rbo1 lc Paraboloidal shell which ha. 
their axes and have ll . Jl " e “ be,s “ r * rl « ld ln direction of 

- proximate fi-iKES ^ ^ , 

.... »,„,r 

v V f -o.fi ^3, 


2 2 

DV V w - 


2 e 

“ab 


jIA & 

+ p o c 


Aj*r 

JxJy 


sir:?",: of 

an eigenvalue problem in which the smalleat .< g 1 proble ® " a ® then replaced by 
shearing stress' at whlc£ The pl^ ITcllZ. —-Ponded to the critical 

•mis work was reported ln full in a published article [ 4 ] . 

K. Lister, Spherical Wave Propagation . 11 - 15 , 12-39. 

used the properties of characteristics to relate the flow at the"“einnti Pr °^ aS 
the end of a given time interval, but did not Involve the of * “ in * “ nd “* 
terlstics as used by Fox and Ralston It h.H ,h , a * r * d °* charac- 

H. Glantz, Aerodynamic meory . 6-13, 7 - 14 , 8-18, 10-13. 

Solutions to the difference equation corresponding to the Laplace equation 
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V U(x,y) = 0 


were obtained In the upper half plane when the following boundary condition* were 
specified: 


11m U(x,y) = 0 

y -«e 

U(x,0) =■ 0 , 


- ?U(x,0) 

jy 


1*1 « i 


lx'* 1. 


Solutions were obtained after reducing the problem to a system of linear algebraic 
equations. A method was also developed for obtaining approximate solutions to the 
more general problem in which the boundar> condition 

-0 . M-n 

was substituted for the second one above. Tlic results arc significant for two- 
dimensional thln-wing aerodynamic theory for tncoopresslble fluid flow. Correspond¬ 
ing three dimensional problems were also investigated. 

A complete write-up of this work is available In the project's third 
technical report [ij. 

L. SeteL, Ap plication of Information Theory to Machine Computations . 11-5, 12-12. 

It was shown In this project that If a measure of gain of Information Is 
associated with the addition of significant digits in a number, then It may pay, 
in the sense of Information carried to the sum or difference of two numbers , to 
retain digits in the summands otherwise discarded. In general, If a sum or differ¬ 
ence is known to lie within an Interval jS rather than 06 ,.^ * 0 ^ , due to the 
addition of significant digits in the summands, then the information gained by using 
the extra digits is defined to be 

loe 

*10 

This definition Is suggested by information theory. This criterion leads one to a 
more meaningful use of the term "significant digit" than Is ordinarily available. 

Similar considerations can be carried to products. 

L. Segel, Quadratic Theory of Bending for a Circular Plate . 12-6. 

A power series solution was obtained for the problem of a circular plate 
under the action of a uniform compression and a uniform bending moment The non¬ 
linear equations and their boundary conditions were non-dimensionalized and then the 
equations were substituted by recurrence relations after expressing the two unknown 
functions In power series. The two unknowns to be solved for were the angle between 
the horizontal and the middle surface and the product of the radius and the horizon¬ 
tal stress resultant. This series approach rather than an asymptotic approach was 
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fr^pSTmS Sp^rhi^dniSr.^r co,n,)utin ‘ *«- 

vert leal deflection a’t ^be 2 ^ ^ ** -‘ tl ° * 

The non-linear equations to be solved were of the form: 

L <P> *Dr P* 

L(x) . - g f 

Where L is the cquldlmenalonal operator 

i.a!.♦ * -*.i 

" ^ ' 

A final report of this work was prepared fsj . 

^ii^SoluUon. of the Di ffusion Equation . 12-34, 14-30, 15-13, 16-25. 
The solution of the heat conduction equation 


jf. £ 


- U U I 




= Q , «t 


. 0< t 


, 0< x<Z 


were substituted by the corresponding difference formulae. Use of the analvtic 
solution was made .0 separate the steady state solution from the transient and t 
determine when the transient became significant crli.n, r transient and to 

the difference approximation were worked out so that romput eY till * n ) S h ilCln8 
mixed and convergence was still Insured nea^ ^ngu^rm« Stab!^? * 


A final report of thia work was prepared [7] 


15-1°0 “le-9?^ ^- cu -^ tlon «* of Structures of Rmlmtlons on rimt. s.,. 

.. .. So “* for *‘la* of K - L - °«vl. were used to compute ho- many non- 

l.omorphlc^, ructure. of relationship, may be defined among member, of a finite 
* el . Th * Whirlwind I computer ... u.ed and of special Interest ... the multi- 
length* 1 ^ arlth " etlC P ro « r — *«re needed to manipulate number, of artltrarv 


14 


15 
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Mure |i reel sols , It there are n object* a dyadic relation a»on« these ob¬ 
jects specified b\ an nxn matrix of Is and Os, a 1 m the 1J place Indicating 
jecis specnieu rplattons are non-Isomorphic if 

that element i bears the relationship to J . Tw 

one cannot be obtained from the other by renumbering tue objects. 

The project used much of the basic logic of the Whlri.ind co-puter^ 
operations were conducted In Integers with no round-off and during dl ^“" ““ 
malnders were kept. As a result, tables were prepared for n ranging up [° - .Z* 

exact numbers as high as lifi". These tables displayed the number of reflexive, sym¬ 
metric, ,rreflexive aynmeirlc. am. asymmetric relations. The results were compared 
with asymptotic formulae. 




_ i _ _. 




B. Rankin, The Numi,.•■■Distributi on of Electro ns and Photo ns In Cascade. 13-48, 14-11 
18-11. 


B, Rankin, • . 1tipJe Prediction. 16-7. 

Tie work of N. Wieier on the prediction of single-time series and his 
work on matrix factorisation ..ore combined and extended to bring explicit formulae 
for the piodictlon of any of a sot of n time series based on the behavior of all 
the throughout the past. The series were assumed completely stationary and 

the pre action operator was assumed linear. Formulae for both the prediction op¬ 
erator and tJ. error of prediction were obtained. 

The form of the solution was a series involving the Fourier coefficients 
of the auto-and ^ross-correlations of the time series involved. Significant de 
partures from He form of solution already known for the one series case were 
found. Though all j-iiipu tat Iona 1 difficulties were not investigated, the formulae 
appeared to submit to calculation with reasonable ease for three series. 

The only report given on this subject was considered final fioj . 


Results that had been obtained by the author in abstract stochastic pro¬ 
cess theorv were carried into the formulation of computing the number of electrons 
and photons that emerge Iron a given thickness of matter when a single high energy 
electron or photon enters the matter. The theory used suggested certain simplifying 
physical assumptions and these were found to permit an explicit algebraic expression 
for the quantity desired. Tills latter expression, though In closed fora, contained 
mam terms which were to be added. The latter operations were prograrmed foi com¬ 
putation on the Whirlwind I computer with the aid of the programming algebra men- 
tloned in the following report. 

In the course of the programming certain approximating formulae were ob¬ 
tained for the convolution of many exponential terms and a Monte Carlo technique 
was developed for selecting distinct, connected paths from an array of points. Com¬ 
plete numerical results are yet to be obtained. 

B. Kan kin, The Basic Problem of Numerical Analy sis Expressed in the Language of 
Computing Machines . 15-7, 18-7, 19-7. 

The three reports cited as references for this project consituted a final 
and comprehensive report on the subject. In this project the mathematical ingre¬ 
dients of a computoi were defined in general and a systematic way to select an op¬ 
timum program in the sense of time, for a given problem was investigated. In order 
to avoid impracticality In the use of too general terms, the general problem, once 
it was formulated, was replaced by the form that it might take for a specific com¬ 
puter. The cot.ipuiur, though l»agina:y, could be simulated by any modern one. An al¬ 
gebra for programming on this ideal computer was constructed and it was shown how a 
program once written in thl* algebra could be optimi^d in "a! i ,v.*vn 

The algebra used as illustration In this general problem was found to be 
of rather general aid in programming specific problems irregardless of optimization 
in time. The translation programs from this algebra to the codes of any modern 
computer were discussed. Certain points of basic interest in the logic of computer 
design were also touched upon. 

The results have been published In abstract form { 9j . 


B. Rankin, Deterministic Processes and Stochastic Process es. 20-7 

It was shown how the simple measurable, metrically transitive, and mixing 
transformation defined by E. Hopf could be used to construct random sequences, 
simple values of sequences of random variables, and, more generally, sample curves 
of stochastic processes. The transformation, used as an example in ergodic theory, 
cairips the unit square into itself in such a way that the successive images of a 
point transformed repeatedly are found within a specified measurable subset of the 
square with a relative frequency that is proportional to the area of the subset. 

Tills is the metrically transitive property and associates relative frequency or 
time avo-ages with probability. The same transformation deforms a specified sub¬ 
set asymptotically in such a way that its overlap with any other specified subset 
is proportional to the product of the areas. This is the mixing property and is 
intimately related to statistical independence. 

It was shown how a transformation with the above properties, essentially 
a deterministic mechanism, can be used in the basic foundations of probability 
theory, thus leading to the conclusion,the probability theory is adequately inter 
preted as an asymptotic theory of deterministic processes. 

Because of the fact that the Hopf transformation can be simply defined 
in terms of snift operations on the binary representation of numbers, an implication 
is immediately gained for the construction of random numbers within a digital com¬ 
puter. This latter point was only suggested and is not yet fully developed. 

The single report given was final £llj . 

R eferences to Final Reports , Technical Reports and Published Papers . 

£lj Phyllis Ann Fox, On the Use of Coordinate PerturbatIons in the Solution of 
Physical Problems, Technical Report No. 1 , Proje ct for Machine Methods of 
Computation and Numerical Analysis , M.I.T., Novsm *er (1953). 

(^2J Phyllis Ann Fox, Perturbation Theory of Wave Propagation Based on the 

Method of Characteristics, Jour. Math . and Phys . , 34, (1955) pp. 133-151. 

[3j Anthony Ralston, On Spherical Wave Propagation, Technic1 Repor t No. 2, 
P roject for Machine Methods of Computation ami Numerical Analysis, M.I.T., 

June (1954). 
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Ltj Anthony Ralston, On the Problem of Buckling of a Hyperbolic Paraboloidal 

Shell Loaded by Its Own Weight, Jour. Math, and Phya. , 35 (1956) pp. 53-59. 

[sj Herbert Glantz and Eric Helasner, On the Finite Sum Equations for Boundary 
Value Problems of Partial Difference Equations, Technical Report Ho. 3 , 
Project for Ma chine Methods of Computation and Numerical Analysis , M.I.T., 
October (1954). 

[6] Lee Segel, Quadratic Theory of Bending for a Circular Plate, ^jarterly 
Progress Report No, 1 2, Machine Methods of Computation and Numerical 
Analysis , M.I.T., June (1954) pp. 6-11. 

[7j Philip L. Phipps, Machine Solution of the Diffusion Equation, Quarterly 
Progress Report No. 16 , Machine Method s of Compu tation and Numerical 
Analysis , M.I.T., June (1955), pp. 25-29. 

(sj M. Douglas Mcllroy, Calculation of Numbers of Structures of Relations on 

Finite Sets, Quarterly Progress Report No. 17 , Machine Methods of Computation 
and Numerical Analysis , M.I.T., September (1955), pp. 14-24. 

[ 9 } Bayard Rankin, An Algebra for Minimizing Machine Performance Time Within a 
Class of Algebraically Equivalent Programs, Bull. Amer. Math. Soc., 62, 

July (1956). 

Bayard Rankin, Multiple Prediction, Quarterly Progress Report No. 16, 

Macblne Metho ds of Computation and Numeric.’ Analysis, M.I.T.. June (7955) 

PP. 7-9. 

(llj Bayard Rankin, Deterministic Processes and Stochastic Processes, " 

Quarterly Progress Report No, 20 , Machine Methods of Computation and Numerical 
Analysis . M.I.T., June (1956) pp. 7-8. 


Part II 

Project Whirlwind 
1. REVIEW AND PROBLEM INDEX 


- 


or Whirlwind comuuter . .. - ! 01 :,i:8 • 8^, hour * ot the 818.07 hours 

(Sf.EC) Group. T he 6 1* 07 incited S % lentl,lc *" d peering Computation 
due .0 computer malfunction! b °“ r ' computer time, 2.3% ... down time 

used for terminal equipment teati ° 7j ' 15 * >our ® ot th * allocated time were 

^ - ---onTr^eT^^^^ - -rlou. 

lions. The^result*°of°29 oTth^ 1 " 8 ^ 1 ’ r rt8r C ° Ver — 16 fleld8 «* *PPUc- 
theses In th^e 29 LobleL V ‘ h * V<? bee " ° r * U1 be ‘"^uded In ac.d. ic 

and 2 Bachelor's Twenty 'tVve 5%'! 2r * " eprMentod 21 doctoral theses. 8 Master s, 

..... .. ..- 

=SSSSrSSF 


Bayard Rankin 
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PROBLEM INDEX 


Aaruatullul 

tallRMriof 


Aerophystce 
■eeeerrb Group 


Chemloal 
Engineering 


lUnrictl 
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Knar ( y 
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Litxi ratory 
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laboratory 


Ibiclear 
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Office of 
Iteiiettcal 

Services 


•ervomecbsalai 

Laboratory 


Spect rosrovv 

laboratory 


Transient response it aircraft structures to aerodynamic beat lac 
Nor Isontal stabiliser analysts 
Paraaetrlc study of coupling and Heapla* 

Hlasi response of rotor blade* 

Trajectory calculation* for a rocset tlur log powered flight 
Lai ran toe of stability derivative* frn flight teat data 

Critical seas calculations for cylindrical geometry 


Nespuna* of a aultt story fraee building under Jynaalr loading 
usepuns* of a Slagle story concrete building tu dynast. loading 
Uynaalc response of shear nails 
Unsign of spherical shell angeenia 

weapons* of a reinforced concrete roof systen to dyaanic loading 

1‘roductlon for trennportetton probien 
Snitching circultn 
Calculation of pi inn uuabeis 

Op11nian11on of alternator control ayatee 


Fourier ayntheaia for crystal atnicturaa 
Fourier ayntheaia for crystal structure 


Stokes particle velocities 


eigenvalue probien for propagation of I.M eavee 
Tropospheric propagation 
Krror analyela 
Prediction analysis 

Arnoapberlc propegatioe of radio savea 
Covsrag* analyela 

rhamodrnaair and dyaanic affects of enter Injection Into high 
tenpersture high velocity gee atreans 
Leninar boundary layar of a steady, conpreaatbl* floe tn the 
entrance region of a tube 
Rolling bearings 
Oroslh of fatigue cracha 

Aaynpi.it I c integration of aquations concerning Torroidni shells 

Spectrsl analysis of atmospheric data 
Statistical end dynamic methods tn forecasting 
■eattier prediction 

Compilation of variances and covariances 

Cal cut at tone for the HIT reactor 
Determination of critical mass 
Flu* leveling In hunogeneoua reactor pert I 

Flua leveling in n<im.ig*n*uu» reactor part II 

Finding eigenvalues of an aryavlric net r Is 
Tenpersture distribution In a bees 


Coulomb va-e function* 

Datsminatlon of pnas* shifts trim experimental croea sections 
Knergy bends in graphite 

Rxchaoge integrals beteeen real slater orbitals 
Neutron-deuterun ecntterlng 
fbeory of neutron roactione 

Al-w sa applied to fao* and body centered iron 
Knergy levels of diatonic hydrides 
evaluation <t teo center molecular Integrals 
Analysis of air shover data 

Iftilllple scattering of eaves from a spatial srray of aplterleal 

acattarara 

Rnsrgy levels of diatonic hydrides LIN 
Atonic eave fund tons 

Polarisabiltty effects In atone and nolaculea 
•Mr* and lapure potaaalun chloride crystal 


taste reduction 
Flight simulation 


Crmple* sped run ennlyala 


Ultrseontc delay tinea 
•Wl 1103 translation prngran 


Minutes of Supervisor or 
•■I Tin> Prog rammer 

SO.I L. Schell 

170N.S C. Crier lone 

178 H R. Patmore 

lil.l K Foes 

Ills * J Prtgge 

/*o.a L. L. Mass..is 


R. 0. Gray 
B. Landry 
Shut Ho 

I Traum 

B. Landry 


J. Dennis 
C Roth 
N Cohen 


■ J. Buergei 
R. MlIsekt 


H R. Delgbt 

P Meson 

L Pataraoa 
L. Peterson 

• Mason 

L. Peterson 


T. T. Toong 

A. She shety 
J. B Saleh 

«. H Crandall 

B. Saltsnan 
l Knlley 
J. Austin 

D Gilman 

T. Cantwell 
J, Barnett 
ft. Kennedy 
k. Hoover 

J Hoseasn 

P. M. Versub 


A, Tubls 
F. Eppllng 
F. J. Corbeto 
P Merrymen 
L. Sartor 1 
N Peshbecb 
J. Pood 

O. Koster 

H. AgheJen lan 
f. Scberb 

P. M. Mors# 

A. Kero 

ft Neebet 

L C. Allen 
L, P. Noel and 

D. T. Rosa 
H Connelly 


0. Arenberg 
J. Pmnkovlch 


•Utunmet 1 cal Problem 

1 ■ *®eiria algebra and equations 

*<x>t of determinantal equation 
Linear equations 
Eigenvalues 
Eigenvalues 

urtl.onormnl last ton. eigenvalues 

Inversion 

Linear aquatluna 

Diagonal last Ion, inversion 

Inversion 

Arithmetic opera*. 1 uaa 
Inversion, eigenvalues 
Eigenvalues factoring polynomials 
Solution of determinants 
Arltnmettc relation# 

Boolean algebra. Hoolear. matrices 
Arithmetic operations 

2 Ordinal) differential aquations 

leven non-lines- first order 

S stem 

•quBlloMs. 

#sv* euqatlon 

Different sets of differential equations 

Second order equstiona 

Simultaneous differential equal tuna 

S*t. of ..•.Him >«»<« >ita*i dlff.r.nll.l 

lair of aimiil tannoua second order 

•Aet of differential equstiona 

Second order 

Set of differential aquations 

3. Partial differential equations 
Second order parabolic 
Second order, variable coefficients 

A. Integration 

Stationary point of a variational 

overlap integrals 

Overlap integrals 

Integrations 

Prssmel integral 

Integrations 
IntegratIons 
Integrat Ion 
Integration 
Integrations 
Integration 

3 Transcendental equations 
Curve fitting 
Nan-linear aquations 
Non 1 loser equations 
Curve fitting 
Non linear aquations 

6. Camples algebra 

Complex roots and function evaluation 

7. Data reduction 

data reduction 

8. Fourier series 

Fourier nyntiisnls 

Sunning serins 

Fourier synthesis 
Fourier summation 

9. Statist ten 

Evaluation of covariance and variance 

10. Integral equations 
Volten* type 


Iteration 

Croul 

Diagonal laat Ion 
Diagonal 1sat Ion 

Scheldt process, diagonal last Ion 

Crout 

Croat 


.»u resslvs sppfos lea I Iona, trout 

Crout 

Direct 

Crout, iteration 
Frame Ultcncom 
Cramer s rule 
Direct 

M> riuakav <*itn» rawed lun 


Fourth order Kutta Gill 
0111 

•rdar Difference equatlona 

Milne predictor-corrector formula 

Step t*v atep 

Gill a method 

Finite difference 

Beckeard difference equal tuns 

Finite difference 

Kutta-ClII 

Finite difference 

Kutta-Gill 


Finite difference 
Finite difference, trapesoldal 
integration 


Simpaon a rule 
Evaluation uf analytic forme 
Evaluation of analytic forma 
Simpaon a rule 

Conversion power aerlna with a 
complex arguaent 
Barnet and Cmilsoti expansion 
Gausa quadrature 
Finite difference 
Irapexoldal rule 
Series substitution 
Finite difference 


Least squares 

Steepest descent 
Iteration 

Least squares 

Staapaat descant 


Polynoniel fitting, ate 


Direct evaluation 
Dlracl evaluation 
Direct evaluation 
Direct evaluation 


T*ble 2-1 Current Problem. Arranged According to Field of Application 


Table 2-11 Current Problem. Arranged According to the Mathematics Involved 


20 


21 




APPROVED FOR PUBLIC RELEASE. CASE 06-1104 


2 WHIRLWIND CODING AND APPLICATIONS 


2.1 Introduction 

. ,4 th« v.rlou* prograaaen are preaentod In 

Programs rapor • am u . .—, , y report preaenta the coatjlned efforts 

numerical order In Section . •' ' , rnhlcis undertaken by members of the Machine 

of Die Projects 6345 and 6915, reports on problems “ n “ erta ~ > 2 2 of Part II to 

Methods of Computation (MMC) Group have been omitted fro. " Group problem, 

avoid duplication of Part I For reference purposes, a Hat of the MMC Group problem. 

appears below. 

Letter, have been added to the problem number, to Indicate whether the 
problem la for academic credit and whether It is sponsored The letter, have the 
following significance. 

A Implies the problem la NOT for academic credit, la DNaponaored 

B Implies the problem IS for academic credit, la UNsponsored. 

C implies the problem la NOT for academic credit, IS sponsored. 

D Implies the problem for academic credit, IS sponsored. 

N implies the problem Is sponsored by the Office of Naval Research. 

L implies the problem is sponsored by Lincoln Laboratory. 

The absence of a letter Indicates that the problem originated within the SfcEC Group. 
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2.2 Progress Reports 

100. COMPREHENSIVE SYSTEM OF SERVICE ROUTINES 


The Comprehensive System of Service Routines refers to a collection of 
utility programs which have been in routine use at the Digital Computer Laboratory 
for the past several years. Included among these are a translator-compiler, a 
post-mortem program and a utility control program which supervises the operation of 
the Whirlwind computer. 


(Airing the past year no significant changes have been made in the vocab¬ 
ulary of the compiler. Some development has taken place, however, in the utility 
control program and the post-mortem program. 

(1) Director Tapes : 

Normal operation of the WW1 computer is initiated by submitting a per¬ 
formance request to the operator. The performance request lists the tapes re¬ 
quired for the run and their operating instructions. A standard set of abbreviations 
has been adopted for use here which conveys all of the required information to the 
operator. The purpose of the director tape Is to replace the computer operator in 
the system. A director tape is prepared by typing the operating Instructions from 
the performance request on a flexowriter. The run can then be executed automati¬ 
cally by reading into the computer the director tape and the tapes forming the run 
(suitably spliced). 


Director tapes have not been as widely used as was expected. There are 
several reasons for this: 


(a) Learning to make use of director tapes places an additional burden on 
the programmer. Since most runs are fairly simple, most programmers have not felt 
any need to make use of the system. 

(b) The system as presently constructed attempts to duplicate exactly the 
manual operation of the computer and requires that paper tapes forming a run be 
spliced together. This has tended to make the system rather unpopular with the 
tape room. 


Director tapes are being used, however, by a small number of programmers 
whose runs are complicated or out of the ordinary since they guarantee that the 
set of operating instructions required for their run will be performed correctly 

(2) The Post-Mo r tem Program : 

A new generalized post-mortem program is being added but is not yet in 
routine use. 

Request for print-outs will be made using a request tape as before. The 
programmer may request that specified ranges of storage be printed-out in specified 
forms (instructions, floating-point numbers, Integers, etc.). Both range specifi¬ 
cation and printed output may be made in terms of floating addresses (previously all 
information was numerical). 

Post-mortem can be obtained as before by reading in a request tape after 
a run has been executed? however, in the new program the following option is also 
available: the programmer may specify up to 7 different intermediate break points 
in his program s.id associate post-mortem requests with each of these. If any of 
the break points are reached during the operation of the program the required post¬ 
mortem will be given and operation will continue. A final break point (with 
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if a program does not terminate at the 

associated requests) can also be apecl . be specified. Requests 

final break point an in case of trouble post-mortem.can P 

using break points must be read in before a program is operated. 

F. Mel wig 

Digital Computer Laboratory 

The CS II conversion routine ha, been modified to accept -muonic c ‘ * 8 
the binary code UOlll in order to make it easier for I « BoJ1 fiction pre 

struction and the 4192 more register. of '°™ ex . clVy 2 S« finds in s tape to 

vents an anomalous stop in conversion when there conversion 

be converted. Previously, more than 25b flads were necessary , or 256 

poat-mortem print too many flads . This post mortem will now be given for M 
flads. The maximum of 255 flads In an acceptable program has "^ been changed^ 

These modi! icatlons have been made in the paper tape of the program and will b 
transferred to magnetic tape unit 0, at the next opportunity. 

A routine for printing out utility programs in parallel columns of octal 
fractions and Whirlwind instructions, with core memory addresses and drum addresses 
has been developed. It la used to provide a format suited to annotation and then 
duplication, so that up-to-date annotated copies of the basic utility programs , e.g., 
the Utility Control Program, the Generalized Post-Mortem Program, may be available. 

The output of the program consists of punched paper tapes which can be run through 
a Flexowrlter for prlnt-out. 

M. Weinstein 

Digital Computer Laboratory 

120 B.N. THE AEROTHERMOPRESSOR 

During the past quarter, theoretical computations were made to check ex¬ 
perimental results. It is expected that the problem will be completed by the end 
of 1956. 

A. J. Erickson 

Mechanical Engineering 

122 N. COULOMB WAVE FUNCTIONS 

The computation program for Coulomb Functions has been terminated tempo¬ 
rarily since the present writer has left the machine computation group. The program 
described in the June 15, 1956 Quarterly Progress Report of the computation group 
has been used by the author this summer while at the Los Alamos, New Mexico, Scien¬ 
tific Laboratory. Tables of Fo, Go and their derivatives for = 4.0, 4.2, 4.4. 

10.0 and P - 0, 0,2,. 40 were computed on the l.D.M. 704 with an accuracy 

between 5 to 7 significant digits. The computations could easily be extended to 
include the range 0 • 7 £ 4 and 10**1* 15 . 

A. Tubis 

Physics Department 


126 C. DATA REDUCTION 


WHIRLWIND CODING AND APPLICATIONS 

Servant.**banisms Laboratory. The over-all problem is composed of many component sec¬ 

tions which nave been developed separately and are now being combined into complete 
prototype programs, inscriptions of the various component sections have appeared in 
past quarterly reports. After tho development and testing of the prototype Whirl¬ 
wind programs is completed, the programs will be re-coded for other, commercially 
available, 1 urge scale computers, (probably the EitA 1103, IBM 701, and IBM /04 
computers), for use by interested agencies for actual data reduction at other loca- 
1 ions The programs are currently being developed by Douglas T. Ross, David F. 
McAvinn, Waiter E. Weissblum, Benson U. Schoff, and Dorothy A. Hamilton, Servo¬ 
mechanisms Laboratory staff members with the assistance of John F. Walsh. This 
w rk is sponsored by the Air Force Weapons Guidance Laboratory through Division of 
Sponsored Research Project 7138. 

Because of its nature, the problem requires extreme automation and effi¬ 
ciency in the actual running of the program, plus the presence of human operators 
in the computation loop for the purpose of decision-making and program modification. 
Thus, while the computer is running, extensive use is made of output oscilloscopes, 
so that the computer can communicate with the human and manual intervention regis¬ 
ters and so that, tho human can communicate with the computer in terms of broad 
Ideas. The computer program thon translates these ideas into the detailed steps ne¬ 
cessary for program modification in confirmation of the human operator’s decision. 
The program which does this translation and modification is called the Manual 
Intervention Program (MIV). The most recent version of the prototype data- 
reduction program is culled tho Basic Evaluation Program. 

During the past quarter, work has continued on the Evaluation Program 
und the MIV System. The evaluation program is being greatly expanded to handle 
considerably more complicated data reduction. Several more test runs have been 
made and compared against previously computed values. A special purpose post¬ 
mortem program written some time ago for this problem has been expanded and incor¬ 
porated into the manual intervention system using intervention switches as well as 
the direct input typewriter keyboard. An elaborate program for controlling the 
labelling and format of both the scope plot and the post-mortem programs is now 
being revised. Thi6 program will later be expanded to Include facilities for 
giving automatically edited logs of manual intervention operations as well as auto¬ 
matic simulation of manual actions for automatic playback. 

As mentioned in previous quarterly reports, programs which are being 
developed on Whirl wind Computer are prototype programs which are to serve as models 
fur production Univac Scientific 1103 computer programs. The 1103 Input Transla¬ 
tion Program of Problem 256 has been jointly sponsored by Problem 126 for the 
purpose of aiding the tranjeription of 1103 computer programs. IXiring the past 
quarter, extensive use has been made of this facility and many programs have been 
translated and run on the 1103 computer at the Air Force Armament Center, Eglin Air 
Force Base, Florida. The details of the procedure were worked out in a visit to 
the Armament Center and at present the debugging of programs is being done entire¬ 
ly oy mail and telephone. It is expected that the majority of this work will be 
completed during the next quarter after which Problem 126 will return to its basic 
function of research prog ramming on the application of digital computer techniques 
to selected problems of tho Air Force Weapons Guidance Laboratory. 

D. T. Ross 

Servomechanisms Laboratory 


Problem 126 is n very large da la-redact ion program for use in the 
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162 N. NUCLEAR SCATTERING PHASE SHIFTS 

During this quarter, solutions and phase shift combinations have been 
found at the finer mesh points necessary for prediction purposes. A theory of be¬ 
havior has been formulated so that future work will consist of finding phase shift 
combinations at higher energies (and a few at lower) that are consistent with this 
behavior. 

Programmers working on this problem are E. Mack and E. Compbell. 

E. Mack 

Nuclear Science Laboratory 

172 B,N. ENERGY BANDS IN GRAPHITE 

The energy bands of two-dimensional graphite have been calculated. These 
results, which were reported in essential detail in the previous progress report, 
formed the basis of a Ph.D. thesis, (Physics Dept., Sept. 1956). 

In view of the two-dimensional results, a survey of three-dimensional 
graphite energy bands has become highly desirable. The major obstacle to this 
extension of the calculation is the writing of a new matrix element generation 
program. Current plans aro to construct such a program. 

F. J. Corbato" 

Physics Department 

193 L. EIGENVALUE PROBLEM FOR PROPAGATION OF E. M. WAVES 

Mode sum calculations at 410 Me and 3000 Me, using the bilinear model, 
have been extended to higher mode numbers. 

Eigenvalues and eigenfunctions have been obtained for an inverse-square 
model at 50 Me, with 500 and 30 foot antenna heights, using both power series and 
asymptotic series. Mode sums were also computed. 

It is planned to continue these calculations until convergence of the 
mode sums is obtained. 

Programmers working on this problem are H.B. Dwight, R.M. Ring and W C 

Mason. 

ft. M. Ring 
Lincoln Laboratory 

199 N ‘ !^IONOF°A N TOB^ LAYER 0F A STKADY ’ COI ®' liESSIBLB FL0 * IN THE ENTRANCE 

A program is being prepared for the solution of the third set of differ¬ 
ential equations for the case where temperature dependence of the vl^cLity a^d 
thermal conductivity of the compressible fluid is taken into consideration. 

Solutions were obtained for the ease »he-e , 

speed is heated at the tube wall. flowing at supersonic 


T. Y. Toong 

Mechanical Engineering 
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...»... 

be modified accordingly to the a, . 8 11 as a ba9lc program which will 

^r - - - 

Programmers working on this problem are R.G. Gray. Shui Ho and E.A. Lawlor. 

H. G. Gray 

Civil and Sanitary Engineering Dept. 
204 N. EXCHANGE INTEG.LALS BETWEEN REAL SLATEH ORBITALS 

parameter range will be notablv ..!«• a / -It is expected that the useful 

Will appear u the “ext progress report with V7 «* iteration 

and subsequent operation. ‘rotary of the results of tostlng 

P. Merryman 

Laboratory of Molecular 
Structure and Spectra 
Department of Physics 
University of Chicago 

216 C. ULTRASONIC DELAY LINES 

time to tl Ie? Sr ‘ Um " ine b ° 0 " COnpleted and d * 9 Kr» 9 will be worked on from 

Programmers working on this problem are Richard Bishop and John Ackley. 

R. Bishop 

Arenborg Ultrasonic Laboratory, Inc. 

225 H,N. NEUTRON-DEUTERON SCATTERING 

author's ‘ "" calculatl ”"» "^ch formed the basis of the 

author s Ph.D. thesis have been extended through program testing The res.il.. 

these revised calculations should be completed ^ affir^ d^ng llT^rL. 

Programmers working on this problem are L. Sartori and H. Paul. 

L. Sartori 
Physics Department 

236 C. TRANSIENT RESPONSE OF AIRCRAFT STRUCTURES TO AERODYNAMIC HEATING 

an" y %r:n h e s°:: : r ?:oL7r„:r erature di,tribution probi * ma ^ ' 

blem will * limited"'^: “ PPe ‘ r “ fUtlJr ' ““ ° f th * WWI Co - pu ‘« r °» this pro- 
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a) production runs of existing programs, 

L) determination of characteristic values and 
c) Inversion and multiplication of matrices. 


vectors 


fur matrix equations, 


L. A. Schmlt 

Aero-elastlc and Structures 
Research Laboratory 


245 N. THEORY OK NEUTRON REACTIONS 

The calculations of the cross sections for c=l.65 and / =.08 wore com- 
pieted, For these same parameter, the angular distributions were done for potential 
wells of depth 38 Mev. and 12 Mev. These results were presented at the Amsterdam 
meeting on Nuclear Physics In July bv Professor Welsskopf. Since then the angular 
distributions for 4.1 Mdv. have been calculated. At the moment small test runs are 
being made with various values of c and J in hopes of getting a better fit to the 
experimental data for the P wave maximum around x = 9.0. If this is possible another 
case will be done for tho values of c and ¥ which give the best fit there and at zero 
energy. 


WHIRLWIND CODING AND APPLICATIONS 

Programmers working on this problem are N Hobbs and E. S. Crlsclone 

E. S. Crlsclone 
Aero-elastic Laboratory 


260 N. ENERGY LEVELS OP DIATOMIC HYDRIDES 

Tho dipole moment of OH was calculated for three values of the inter* 
nuclear distance An electron distribution analysis was carried out using the re¬ 
sults of the SCF-MO calculation described earlier. 

The problem has been terminated. The results are being written up for 
publication 


A. J. Freeman 
Solid State and Molecular 
Theory Group 

261 C. FOURIER SYNTHESIS FOR CRYSTAL STRUCTURES 


E. Campbell 
Nuclear Science 

253 N. APW AS APPLIED TO FACE- AND BODY-CENTERED IRON 

Little machine work has been done on the problem during the last quarter. 
At present the particular electronic states at the k point (0,0,0) are being inves- 
tigated. 


The crystal structure investigations of pectollte, dlglycine hydrochloride. 
and dlglycine hydrobromide, mentioned in earlier reports, are complete and have been 
submitted for publication. 

The refinement of the structure of wollastonlte has been successfully 
continued and the structure of rhodonite is under investigation. 

Programmers working on this problem are M. J. Buerger and T. Hahn. 


Programmers working on this problem are J. Wood and M. Saffron. 

J. H. Wood 

Solid State and Molecular Theory Group 

256 C. WWI-1103 TRANSLATION PROGRAM 

A translation system for the ERA 1103 computer has been developed under 
a contract with the Servomechanisms Laboratory at M.I.T. With this system a pro¬ 
grammer can write programs for the 1103 Computer (in a specified language closely 
resembling that used in the comprehensive system compiler) and have them trans¬ 
lated by WW1 into binary tapes which are acceptable to the 1103. The program is 
working satisfactorily at present. 

F. Helwig 

Digital Computer Laboratory 

257 C. HORIZONTAL STABILIZER ANALYSIS 

In this quarter the complete program which analyzes the response of a 
present day fighter to a gust type disturbance has been tested. The analysis em¬ 
ploys numerical integration techniques to solve twelve non-linear differential 
equations of motion. At tho present timo a study Is under way to assess the im¬ 
portance of some of the parameters associated with the problem. Within the next 
quarter it is anticipated that many production runs will be made in an effort to 
complete the problem. 


M, J. Buerger 
Geology Department 

262 N EVALUATION OF TWO-CENTER MOLECULAR INTEGRALS 

The values of hybrid Integrals evaluated by Dr. Corbato's routine were 
checked against Lh© values of the same integrals obtained elsewhere and they were 
found to agree to four decimal places. This accuracy is not enough for the spe* 
clflc problem which is being calculated and therefore, attempts are being made to 
evaluate these integrals by P. Merrymans routine. 

H. Aghajanlan 
Solid State and Molecular 
Theory Group 

264 C. OPTIMIZATION OF ALTERNATOR CONTROL SYSTEM 

A program has been written to perform simple optimization with respect 
to the weight of an aircraft alternator. The optimization problem was simply 
phrased for the first computer runs in order to provide information regarding the 
operation of the program when the functions involved are highly nonlinear. The 
problem was phrased as follows: minimize the weight of the alternator holding 
the full load, zero power factor losses constant, by varying the magnetic flux den 
sity in the air gap and the armature stack length. This simple, 2-variable, one- 
equality constraint problem was run successfully for several different terminal 


28 


29 









APPROVED FOR PUBLIC RELEASE. CASE 06-1104. 




WHIRLWIND CODING AND APPLICATIONS 

condition*. Tlic results o I Ihi-se opt U.li»t Ion runs ln<Jar equality 

vcloped for minimising * nonlinear function 1 P - operate successfully, 

constraint* by the method of steepest descent can be mao 

A complete description of the techniques used cm be found 

183 Deslan of Systems Using Digital Computers , (to be published). 

technique is applied to a practical design problem. These difficulties are b lng 
investigated In order to minimize them and to extend the class of problems which 
can be solved by the technique. 
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274 N MULTIPLE SCATTERING OK WAVES FROM A SPATIAL ARHAY OK SPHERICAL 
SCATTER BUS 

All values needed tri compute the nine Hf expressions for seven directions 
have been obtained. 

A profcraa is now being tested to fit polynomials to the My s by the me¬ 
thod of least squares. 

This is being done to facilitate Interpolating for the Hf's. 


It is anticipated that an optimisation of the aircraft alternator sub¬ 
ject to several constraints will be run In the ensuing quarter. The r«»dUs of 
such a run should provide practical information regarding the physical and per¬ 
formance characteristics of the optimum machine. 


Programmers working on this problem are J. B. Dennis and R. R. 


Brown . 


R. R. Brown 

Energy Conversion Laboratory 


266 A. CALCULATION KOR THE MIT REACTOR 

Since the design is now firm for the M.I.T. Reactor, the core was cal¬ 
culated using these final properties. A new type of core, with only 20% U-235 in 
the fuel, was also calculated. 


M Kurakashlan 

Nuclear Science Laboratory 


278 N ENERGY LEVELS OF DIATOMIC HYDRIDES Llll 

IXirlng this quarter we have extended the calculation of the electronic 
energy of lithium hydride to lnternuclear distances of 7 and 8 atomic units, using 
the non-orthogonal atomic orbitals and the Lbwdln formalism described in previous 
progress reports. The potential energy values which were derived from these re¬ 
sults for the ground stale and the first excited state are consistent with our pre¬ 
vious work and the potential energy curves drawn through all the points show no 
discontinuities. 

A self-consistent field calculation and a subsequent configuration in¬ 
teraction treatment using orthogonal molecular orbitals were carried through for 
the 9 lnternuclear distances at which we have evaluated the electronic energy. 

This second procedure shows remarkable agreement with the non-orthogonal case. 


At the moment no future calculations are planned except for those com¬ 
pleting the above sets which are 90% complete now. 

Programmers working on this problem are F. Bradley, and J. O'Connor. 

T. Cantwell 
Chemical Engineering 

270 B. CRITICAL MASS CALCULATIONS FOR CYLINDRICAL GEOMETRY 

During the past quarter, all program work was finished. The results 
are currently being drafted Into a thesis: Nuclear Characteristic of Heavy 
Water Moderated , Enriched , Homogenous Reactors , which will appear shortly. Read¬ 
ers arc referred to this report for full details of the programs and results 
achieved. 

James R. Powell 
Chemical Engineering 

273 N. ANALYSIS OF AIR SHOWER DATA 

Hiring the past quarter we have used 7 hours of computer time to ana¬ 
lyse 500 events and to do a Monte Carlo study of the resolution of the experiment 
Wo plan to continue the analysis of data and finish the Monte Carlo calculations. 
A complete description of the program and method of analysis will be written in 
the next month. 


Using the molecular wave functions which have now been evaluated, work 
can be continued in evaluating other molecular quantities such as the dipole mo¬ 
ment. A more detailed description of these recent calculations is given In the 
current (Quarterly Progress Report of the Solid State and Molecular Theory Group. 

Programmers working on this problem are A. M. Karo and A. R. Olson. 

A M. Karo 

Solid State and Molecular 
Theory Group 

288 N. ATOMIC WAVE FUNCTIONS 

Calculations for various atomic functious are being carried out. such 

as: 

1) He conftguiation interaction 

2 ) transition element calculation 

3) Wave functions for an MnO band calculation 

4) four electron systems 

Professor Nesbet's program has been used to do these first three sets of 
SCF calculations: 

1) In sn atomic He where 20 configurations are taken in Interaction 
(Corbato's matrix diagonallzatlon 1* also used here.) 
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F. Scherb 

Cosmic Ray Laboratory 
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2 ) in various configurations of transitions metals In their atomic state, 

such as Fes, (15) Z (25) 2 (2p ) 6 (35 ) 2 <3p ) 6 (4s ) 2 and (15) 2 (25) 2 (2p ) 6 <35 ) 2 (3p ) 6 

(3d ) 7 (4s ) 1 and (15) 2 (29) 2 (2p ) 6 <35 ) 2 (3p) b (3d) 8 . 

3) In calculating M„++ and 0* Ionic functions, which will be used as the 
basis of an energy band calculation for MnO. (This last calculation will be done 
under a new problem number.) 

4) Calculations on four electron systems, L1H, He2, and Be have been com¬ 
pleted. Calculation of the magnetic hyperflne coupling constant of LI has been com¬ 
pleted. Calculation of the binding energy of CH 4 by the one-center method, with d 
and f orbitals, Is In progress. 

During the past quarter, the He work was essentially completed. The 10 
Ingoing parameters were varied In order to produce the best final energy. Numerical 

discrepancies were found In some of the results of our transition element work. IXir- 

ing the next quarter we hope to complete this work. 

The MqO problem is a major calculation, which has Just been started under 
another problem. 

Programmers working on this problem are R. K. Nesbet and R. Watson. 

H. Watson 

R. Nesbet 

Physics Department 

290 N. POLARIZABILITY EFFECTS IN ATOMS AND MOLECULES 

The Whirlwind techniques now being used are the same as those described In 
previous reports. However, during the last quarter a great number of results of phy¬ 
sical Interest have been obtained. These are described In the Solid State and Mole¬ 
cular Theory »iarterly Progress Report for October 15, 1956. " 

L. C. Allen 

Solid State and Molecular 
Theory Group 

293 C. ROLLING BEARINGS 

A new Iterative gradient method for solving non-linear simultaneous equa¬ 
tions was programmed and tested. Although successful when used on a sample set of 
equations, it gave negative results when used on the bearings equations, because of 
the ill-conditioned nature of the latter. No further computer work has been under- 
taken. 


A. J. Shashaty 
Mechanical Engineering 


300 L. TROPOSPHERIC PROPAGATION 


“ " a * f °“ nd de8lrable to introduce a third type expansion Involving ex¬ 
ponentials which would be valid for intermediate values of the argument. A. the pro- 

!n?e™ediat and Y *' U86d ^ ^ a sum of exponential, for 

Intermediate value, and an asymtotlc expansion for large argument. 
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C.uy has riulcr^fun; 10 ^.'^.^:* th * — automati- 

sliuation which reduce, to Z Ordinary^ ^ C ~ ° r ub 
With calculation, made by J. A. Greenwood of H^rd agree, l ° flgUrC8 

It 1 . expected that the problem will be completed In the next quarter. 

W. C. Mason 
Lincoln Laboratory 

30fJ U. SPECTRAL ANALYSIS OF ATMOSPHERIC DATA 


309 B.N. PURE AND IMPURE POTASSIUM CHLORIDE CRYSTAL 


B. Saltzaan 
Meteorology Departaent 


prOVlded lntere *tl n g information which should be applicable to 

- r i:ir 

after the principal calculations discussed above. 

L. P. Howland 
Solid State and Molecular 
Theory Group 

310 C. TRAJECTORY CALCULATIONS FOR A ROCKET DURING POWERED FLIGHT 

The entire program has been completely checked during this quarter and 
-sny trajectories have been computed. At pre.en. a parameter study i. ^elng c ‘r 
rled on, which will be continued during the next quarter. * 

2392 ThU W ° rk 18 8pon8 ° re<1 b > tb « u s - **r Force under Contract AF 33(616)- 

J. S. Prlgge 

Aerophyslcs Research Group 

312 L. ERROR ANALYSIS 


There has been a revision In the method of programming the modification 
descrlb^ in gPR No. 46. This revision was nece.s.ry be^.u.e £ dlff^lUe. 
countered In the use of the Programmed Arithmetic. Production runs are being made 
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with this revision. 


Programmers working on this problem are L. Peterson and E. Hutcheson. 

L. Peterson and E. Hutcheson 
Lincoln Laboratory 

317 C. EXTRACTION OF STABILITY DERIVATIVES FROM FLIGHT TEST DATA 

During the past three months Investigation of Shlnbrot's method for the 
extraction of stability derivatives has been continued, covering such aspects of the 
problem as effects of changing the number of frequencies and of varying the assumed 
precision of the simulated response measurements . Determinations of lateral as 
well as longitudinal derivatives have been made. An investigation has been initia¬ 
ted covering a modification of Shlnbrot's method using differences of cosines as method 
functions as well as the sine functions suggested by Shlnbrot. Also, a new method for 
stability derivative extraction based on applying the principle of maximum likelihood 
to an assumed probability distribution for response measurements and substituting 
statistical estimates for the quantities in the resulting equations for which explicit 
analytical expressions do not exist is being considered. In the coming quarter, it 
is intended to continue the investigations of Shlnbrot's method that are presently 
being carried on and to determine whether either of the new methods offers promise of 
better results. For a more detailed description of the work done and proposed on the 
problem, consult the Progress Reports on the Dyn amic Stability and Control of Aircraft 
*nd Missiles of the Aeorphysics Research Group. 

Programmers working on this problem are T. Carney, M. E. Hoult, L. L. 

Mazzola, M N. Springer and L. E. Wilkie. 

M. N Springer 
Aerophyslcs Research Group 

326 C. PRODUCTION FOR TRANSPORTATION PROBLEM 

A number of large transportation problems have been solved during the 
past quarter u.lng the Whirlwind tran.portatIon routine,. These problems concerned 
optimum shipping from 30 to 60 plants to 260 to 291 customer*. 


minute* ^ C °" PUter tiam r *« ulred »°lve these problem, ranged from 30 to 50 

Programmer* working on thl* problem are J. Dennis and W. Jewell. 

(•See description under Problem 219 In Summary Report. 45 and 46.) 


327 L. PREDICTION ANALYSIS 


J. B. Dennis 
Electrical Engineering 


pned to t™: “r thia probieB are now bein * 

method of prediction The effects <r th . ° D * P°»er-»»rles-approxiBation 

also being tested. '"‘h 0 " ^ predlCtl ° n a ~ 

and checked out. An the8e two ®tudies is largely completed 
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is a Maximum Ltkiuh^rf'c’lJer^on’^th^,‘infi d "°' 46 U now bein * Programmed it 
previously programmed Large sections^,progra^ 3 2 Ind Problem 

form a part of this much -ore el.bor.te 0 anai r y*ir “ nd 327 * re b ' ln « c “" b ‘" ed to 

Programming 1. now In progress. 

•m... « «■*■ L. „ 

L. Peterson 
Lincoln Laboratory 

334 C. PARAMETRIC STUDY OF COUPLING AND IWMPING 

planes at various altitudes^nd^lrspeedB h"* 0 h* 6 " co “ pleted on several alr- 

present re^t." W ° r “ °“ tM ' Pr ° ble “ lB d «»»* nd "“ upon a compete evaluation of the 

K. Wetmore 

Aero-elastlc and Structures 
Research Laboratory 

341 C. STATISTICAL AND DYNAMIC FORECASTING METHODS 

binations of the original L»I? “ P W ‘ S undertake "- Those eight llnea, com 

2 222 =" ~ 

The attempt to incorporate temperature parameters into the t.a,. 
tion equation has been completed. The grid field nnt* i . y predic- 

temperature information necessitated the computation of 28°fun«lon« ( 7 “ “ b 5 d,ly " ean 

surface temperature anomalies at several United States 1 * ° 5 * day mean 

-qutr!iru,T“ 1 ^trl t Clnd nd I the ,,PPllC * tl0n ° f the independent "da t a Vm 
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determine empirically the beat values to use for correspondence to actual atmos- 
ph!r™ motions. In order to extend the number of iteration, in them, expe.lment. 
Without the time factor becoming infeasible. . revision of the original program ... 
undertaken. 

All of the work done on the project ha. been supervl.ed by Profeaaor 
Edward N. Lorenz, Department of Meteorology. 

E. A. Kelley 
Meteorology Department 

343 C. WEATHER PREDICTION 

During the last quarter work has continued on the development of routines 
for work on linear prediction for various classes of weather systems. Because of the 
large amount of data to be analyzed, a conversion routine for input of numerical In¬ 
formation only was written. Use of this routine reduces read-in time a significant 
amoun t. 

Future plans call for further development of linear prediction equations 
Quite recently additional data arrived from the National Weather Records Center, 
Asheville, N.C. to further the investigation. 

Work on this program has been under supervision of Professor James M. 
Austin, Meteorology f'epartment. 


Programmers working on this problem are Herb Brun and E. 


Kelley. 


E. A. Kelley 
Meteorology Department 
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354 D. RESPONSE OF A SINGLE STORY CONCRETE BUILDING TO DYNAMIC LOADING 

The original program for this problem was designed to solve a specific pro- 
blem and the required reaulta were obtained. Our recent effort haa been to generalize 
the progrui so that It would be able to handle any structure subjected to eaplrlcal 
loads, provided that the dynamic model was limited to either one or two degreea of 
freedom. The program la ready for debugging and will be used to obtain additional re- 
aulta. 

B. Landry 

Civil and Sanitary Engineering 

360 C. DYNAMIC RESPONSE OF SHEAR WALLS 

lAirlng the paat quarter, production runs were completed, ualng the original 
program. The reaulta are contained in a Sc.D. Thesis, by A. Flnerman, Courae I, 

1956. In the future a new program will be written to account for baae fixity, the 
previoua reaulta being obtained for two-point baae aupport. This work will be car¬ 
ried on by Mr. Shui Ho, who replacea Mr. Flnerman. 

A. Flnerman 

Civil and Sanitary Engineering 

361 B,N. GROWTH OF FATIGUE CRACKS 

In reference 1, the problem of determining the rate of growth of a fatigue 
crack in a rectangular bar compound of an Ideally plaatlc material was analyzed for 
the caae of toralon in the fully-plaatic range. An integral equation expreaaing the 
growth rate was obtained In finite difference form; this equation was programmed for 
Whirlwind, and, during the last quarter, results were obtained for several values of 
the parameter, Xs. 


All available Zeeman data was used to search for energy Intervals. No 
definite Intervals have been found as yet, but several possible ones have been found 
recently. These will be checked In the netr future, and if necessary more Zeeman data 
will be obtained. It is believed that some definite results will be obtained within 
the next quarter. 

J. Lindner 

Spectroscopy Laboiatory 

350 D. COMPUTATION OF VARIANCES AND COVARIANCES 

The set ot eigenvector* mentioned in the last report have been obtained. 

The determination by matrix multiplication of the amplitude coefficients of these 
vectors is the next major task for machine computation. 

Programmers working on this problem are Bernard Shorr and Donald L. Gilman 

D. L. Gilman 
Meteorology Department 


The analysis mentioned above is valid only for cracks originating from a 
sharp notch. The analysis was redone for cracks originating from notches of radius 
X<> and resulted in an integral equation similar to the one obtained in the sharp- 
notch analysis. The new equation was programmed for Whirlwind and is now in the pro¬ 
cess of being corrected. 

1. Quarterly Progress Report No. 20, Mschine Methods of Computation on Numerical 
Analysis. June 15, 1956, pp. 14-16. 


J. B. Walsh 
Mechanical Engineering 


362 B. FOURIER SYNTHESIS FOR CRYSTAL STRUCTURE 


Crystal structure determination and refinement work of the mineral jame- 
aonlte hue been completed. 

N. Nllzekl 
Geology Department 
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363 A. ASYMPTOTIC INTEGRATION Of EQUATIONS CONCERNING TORROIDAL SHELL 

Asymptotic Integration of a pair of alaultaneous second order ordinary dif¬ 
ferential equations was accoaplIsbed for several values of pressure, using finite dif¬ 
ference sethod. The future plan Is to obtain a series solutlou to the problem, tak¬ 
ing enough terms In the series to obtain satisfactory agreement with the asymptotic 
solution. In this manner a solution will be available for all values of the impor¬ 
tant geometric variable. The results to date were reported at the 9th International 
Congress of Applied Mechanics In Brussels during September and were also reported 
before the 11th Annual Meeting of the Association for Computing Machines at Los 
Angeles, In a paper: Optimum Recurrence Formulas for a Fourth Order Parabolic Par¬ 
tial Differential Equation , (August 23-29, 1956). ~ - 

Programmers working on this problem are N. C. Dahl and S. H. Crandall. 

N. C. Dahl 
S. H. Crandall 
Mechanical Engineering 

364 C. BLAST RESPONSE OF ROTOR BLADES 

Ikirlng the past quarter a program has been written to calculate the final 
" atri u.°f th ® flr,t Ph “* e ° f thi “ P 1,01 * 1 ®"- Thi» program, together with a preliminary 
one which computes the elements of the matrices by the solution of difference equations 
is being used with various parameters to solve for the desired normal modes, shapes and 


365. FINDING EIGENVALUES OF AN ASYMMETRIC MATRIX 


M. Callaghan 

Digital Computer Laboratory 


J. Koseaan 

Digital Computer Laboratory 

367 B. DETERMINATION OF CRITICAL MASS 

termlne at what value of uM^Inrichmlnt'ealh'd”!^* 1 *"? W ° re lnve,tl 8 ated to d ®' 
tical mass of uranium. «nrlch«e„t each dealgn would require it. minimum crl- 

An Engineering and EconLlc'clmparlsln t ir n *.r"* 1 ^ C ’!' POr “ ted into * thesis entitled: 
tlons !- - ° — ?«£tor Designs for Pscksge Power Appllcs - 

neering Department with Instructmns^r Ua ^"in^h ° D *" the NUCl * ,r Engl " 
further Investigation In this field * th * vent thero 18 need toT » 


Trooat. 
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Programmers working on this problem are J .. , „ 

“re J. Barnett, J. Powell and M. 




J. W. Barnett 
Nuclear Engineering 
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erly reports. Durl^t^present 'qulrtej*,^^ h “ S b ‘' en glVen ln P r «'« d *nE quart 
tamed for different physical cohfJgur.t^.^^/^^^-l-^ns .e, r op¬ 
tion. for the Whmwir d P pr!grJ! a a8 e «n U ^“ t o 4 cari at0r l ° ob ' 8in check "olu- 

case problems. These solutions were nuir *° ca "» OUt tho sol utlon of two limiting 
both a check solution for the runs obtained on^he^i™? the 8 tand P olnl at providing 
providing useful information for the U^nTc^. Compuisr.aa well as 

seven decimal d lglts*which 1 ls d excerdm“y go^ 65 ° “ nd occa,lonal >V 

terpretlve ayatema were used on SSSS ST 

Computer ll^afry'« T.'S ,£ £ s ° 
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ATMOSPHERIC propagation of RADIO WAVES 


F. M. Venuh 

Office of Statistical Services 


th. Tchbyhhirint.*™!^* d ‘ tB r ° r «t^hph«rle ..fraction 

«• 2 — *_ 

s I'zzzz rr - -••• - —* 

. 1.0 be atmmpted°in°the*future! 0n * ,r °“ 1601 **" -HI 

Programmer, working on this problem are W. Mason and P. Duffy. 

W. Mason 

Lincoln Laboratory 

372 B. DESIGN OF SPHERICAL SHELL SEGMENTS 

During this quarter the membrane stresses of a spherical shell segment 
-ere computed m. future p.rt will deal with bending atres.es snd dUpl.^t.. 

E. Traum 

Civil Engineering 

373 B. FLUX LEVELLING IN HOMOGENEOUS REACTOR, PART I 

Denartwe , re “ Ult * °* thl8 Problem have been submitted to the Chemical Engineering 
Department in a Master s theala entitled Unifo rm Power Distribution In „ H,.fl«et.. ( i 
H omogeneou a Reactor, by L. N. Hoover and R. # Kennedy. 
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For Problea 373 the theory and programing were eupplled by the under 

signed. 


The problea aolutlon on Whirlwind 1 was very satisfactory Results were 
In agreeaent with physical considerations. However, coaparlson with experimental 
data was not possible since such data are not available. 

I would like to express ay appreciation for the outstanding assistance 
of the staff of the Digital Coaputer Laboratory Both tape rooa personnel and 
coaputer operators did their utaost to alnlmlze ay errors, all of which were direc¬ 
tly ay own fault. 
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377 L. COVERAGE ANALYSIS 


The first half of the prograa, as described In QPR No. 46 has been com- 
pleted; production runs have been Bade. The data-handllog aspects of the problea 
are now being considered, Initial programing having been done by R. Karp. 

Programmers working on this problea are L. Peterson, E. Hutcheson and 

R. Karp. 


L. Peterson and B. Hutcheson 
Lincoln Laboratory 


It has been a privilege and pleasure to work with the coaputer laboratory. 


378 B. RESPONSE OP A REINFORCED CONCRETE ROOF SYSTEM TO DYNAMIC LOADING 


R. W. Kennedy 
Nuclear Engineering 


The object of this problem was to study the dynamic response of the roof 
of a single story reinforced concrete building for various loading functions. 


374 B. FLUX LEVELING IN HOMOGENEOUS REACTOR PART II 

A uniform fission rate throughout a nuclear reactor core will Increase 
the performance and efficiency of the reactor by permitting each region of the core 
to operate near Its temperature limit, which Is determined by heat transfer and 
metallurgical properties 

By employing a suitable reflector with a homogeneous core, It is possible 
to greatly flatten the thermal flux and low energy flux distribution in a spherical 
reactor. For an epithermal reactor, this is equivalent to approaching a uniform 
fission rate in the core. 


The problem was the same as 354 except for the loading portion ,and for 
that reason similar portions were not repeated. 

The loading portion of the program consisted of three types of loading 
curves: triangular, square, and sinusoidal. For each of these types the slab and 
beam response was found for three different loading time durations. The results of 
the response of the system for each type of curve were interpolated to find the peak 
value of thm curve of the same shape and loading time duration which would yield a 
slab deflection In a specified range. 

ftirlng the past quarter, this problem has been completed. 


The results of diffusion theory are applied by means of finite difference 
approximations to 1 energy groups and n spatial regions for a reflected spherical 
reactor with homogeneous properties In the core and reflector respectively. 

The finite difference equations solved are: 


— A . A A _ A 1 

F » k F - F - 1 
dtI n n-1 n 


where k 1 and I* are parameters determined by the material properties of the reactor 

The core radius is varied In successive iterations until the critical con- 
dltlon is achieved. 


L H. Hoover 
Nuclear Engineering 

376 N. FLIGHT SIMULATION 


total of 41 apparently correct runs have been made. Each of 


tnese runs 


represent. 10 a.conds of flight of the aircraft. The numerical data 1 . being pro¬ 
cessed, and If It proves to be accurate and conclusive, the problem Is completed. 


H. M. Bourland 
Servomechanisms Laboratory 


B. Landry 

Civil and Sanitary Englnesrlng 


380 B. SWITCHING CIRCUITS 

Two programs have been written for analysis of switching networks by 
means of Boolean matrices. The first of these programs Is for analysis of two- 
terminal contact networks to be analyzed Is represented by s symmetric Boolean matrlz, 
and the terms above the main diagonal are fed Into the computer In a simple coded 
form. After translation of the code, the computer reduces the size of the matrix by 
removal of successive nodes until a two-by-two matrix remains. The standard sum 
which represents the transmission between the output terminals of the network is ob¬ 
tained from this matrix and printed out In decimal form. An 8-varlable, 42-node 
network can be completely analyzed in less than two minutes, and small networks re¬ 
quire only a few seconds. Tbs second of these programs Is for analysts of multiple- 
output contact networks of 4 to 8 variables and 3 to 40 nodes. Tbs terms on and 
above the main diagonal of the symmetric Boolean matrix which represents tbs con¬ 
tact network are fed into the computer, and the computer squares the matrix success¬ 
fully until an output matrix is obtained. The transmissions between any or nil node 
pairs in the network can then be obtained directly from tbis output matrix. 

The McCluskey-Quine reduction procedure Is currently being programmed to 
determine minimum sum forms of switching functions of 15 or fewer vsrlables with a 
maximum of 512 terms In the standard sum. The part of the program which determines - 
the prime laplicanta from the standard sum has been completed and checked, and the 
remainder of the program should be completed In s few weeks. After the performance 
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WHIRLWIND CODING AND APPLICATIONS 

of this program has bean evaluated! a program will be written for ayntheala of non- 
aerles-parallel switching circuits by meana of Boolean aatrices. 

Charles Roth 
Electrical Engineering 

382 B. CALCULATION OF PRIME NUMBERS 
Problem Description 

The feasibility of determining prime numbers on the Whirlwind computer, 
the minimization of the time Involved and the limiting factors are to be determined. 

The sieve of Eratosthenes lists all the numbers to be tested and, then, strikes out 

those divisible by 2, those divisible by 3, those divisible by 5. those divisible 

by pi, where is prime and less than Y, the greatest number to be tested. This 
method will be modified in the following ways in order to make the problem applicable 
to the computer: 

1. Only primes equal to or less than the square root of Y need be considered 
for eliminating composite numbers 

2. X, the least number to be tested, may be reduced mod p, for each elimi¬ 
nation prime, and a least negative residue can be obtained in each case. Since 
Xs^rnod implies that X+lwJr^+fJmod p^, the negative residues can be used for the 

elimination of composite numbers in lieu of the corresponding primes. If at any 
time X+a=0 mod p^ we know that X+a is divisible by Pj^ and hence composite. This 
process of adding one to each negative residue and seeing if it becomes zero may be 
continued until all numbers between X and Y are tested. 

Two variations of the above procedure will be considered to determine 
which will take less time. The first variation checks every residue on each number 
to be tested. The second goes onto the next number to be tested as soon as one 
prime divisor has been determined. Although the second method seems shorter, there 
is more time involved in resetting the residues to congruent negative values. 

Limiting factors such as number of memory storage registers, number of 
bits that can be stored in a single register and economical time will be considered 
in anticipation of further attempts at determining prime numbers on the Whirlwind 
computer. 

Progress During Quarter 

I have divided my problem into four parts, each of which operates inde¬ 
pendently of the others. The first part computes the smallest (in absolute value) 
negative remainders of the fraction, X/p *, for each prime, p I have already ob¬ 
tained a program that will do this for l X's that are small enough to fit in a 
single register, but when this routine is altered slightly to accommodate two regis¬ 
ter numbers it falls. At present I am working on a program that will allow me to 
check the formation of the residues at each step. 

Once the residues have been formed, the smallest prime divisor greater than 
1, of each number n between X and Y is recorded. I have developed a program that 
will do this if given the proper negative residues. I also have a second program 
that can do the same operation, and I intend to determine which of the two is more 
economical. 


WHIRLWIND CODING AND APPLICATIONS 

.. r Th ! 1 thln, J!* rt 0t th * proble " *• recording the least prime divisors in a 
suitable fashion. The outputs that I am using at the present time would not be prac 
tleal for the computation of, say, a million numbers, although it is quite satisfac¬ 
tory for the few hundred numbers that I am using to troubleshoot the program. I have 
not yet begun to work on an output that can handle the large amount of numbers that 
1 will get ultimately. 

The final part, which is dispensable and which will bo completed only if 1 
have time, is the making of a program that will find and store all primes, p s 
This part is not absolutely necessary since the primes could be found in tables and 
recorded on the actual input program. The recording of these primes is somewhat 
laborious and this last part will facilitate further attempts to determine prime 
numbers on the Whirlwind computer. 

H. Cohen 

Electrical Engineering 

383 C. STOKES PARTICLE VELOCITIES 
Problem Description 

The purpose of this problem is the solution of G. B. Stokes' equations 
for the motion of water particles in oscillatory waves of finite amplitude. The 
entire program is being carried out under Contract DIC 7116 sponsored by the 
Humble Oil and Refining Co. of Houston, Texas, and entitled. Wave Forces on 
Offshore Structures. 

The problem is being solved in an effort to compare the theoretical velo¬ 
cities with experimentally measured velocities. The comparison will aid in evalua¬ 
ting the rationality of various relations derived from Newton's Laws for determining 
inertial forces on objects submerged in oscillatory saves. 

The output of Whirlwind Is to be Flexowrlter tape for print out of the 
horizontal and vertical components of particle velocity versus wave phase angle 
The velocities are to be computed for every 15 of phase at three depths for each 
of 12 waves. 

Progress During Quarter 

Tapes have been prepared and run to check the solution. The problem is 
now ready for production runs as soon as experimental parameters are known. It ia 
expected that these parameters will be known in a month and all work will be com¬ 
pleted by November 31, 1956. 

Programmers working on this problem are T. Marlow and J. Hougley. 

T. Marlow 

Hydrodynamics Laboratory 
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APPENDIX 


X. SYSTEMS ENGINEERING 

WW I_K KLLABILITY 

The following is the WWI Coaputer Reliability for the paat quarter: 


Total Coaputer Operating Tlae in Hour* 2572.2 

Total Tlae Lost in Hours 32 9 

Percentage Operating Tiae Usable 98.8 

Average Uninterrupted Operating Tlae 

Between Failure Incidents in Hours 22.3 

Total Nuaber of Failure Incidents 114.0 

Failure Incidents per 24 Hour Day 1.06 

Average Lost Tiae Per Incident In Minutes 17.3 

Average Preventative Maintenance Tlae Per Day In Hours 2.2 

HIGH-SPEED MEMORY 


The storage facilities available at WWI Include aagnetic tape units, aag 
netlc druas, and magnetic core aeaory. The access tiae to both the tape and drus 
storage is long in comparison to that required of the core aeaory. WViI la frequently 
scheduled to provide solutions to many real tiae problems. To meet these demands, 
the high-speed core aeaory facilities have been expanded by a factor of three. The 
addition of 4096 registers raises the high speed core memory complement to 6144 re¬ 
gisters (each register is composed of 17 bits). The WWI coaputer is designed to em¬ 
ploy a 16 bit work in which 5 bits form the instruction and 11 bits address the de 
sired register. Therefore, one-third of the registers (2048) can be actively en¬ 
gaged in problem solving, while the remaining two-thirds are passive. However, 
through use of a new coaputer instruction, change core memory fields (cf), it is pos 
sible to activate or deactivate any of the six core memory fields (a field is com¬ 
posed of 1024 registers) in 20 microseconds. The cf instruction provides several 
flexible features in having coaputer pulses interrogate the address bits for the cf 
instruction to determine if any of the following actions are to be performed: 

1 . program pass control to the succeeding register or to any other 
register; 

2 . replace a single field, both fields, or leave the field assignments 
unchanged; 

3. assign any of the six fields to appear as fixed addressed registers 

0 - 1023 or 1024 2047. 

To further assist WWI users in solving real tlae problems, the core memory access 
time has been reduced from 8 to 7 microseconds. 


2. VISITORS 


„ < Tours the Whirlwind I installation include a showing of the film, 

■taking I-lectrons Count", a computer demonstration, and an Informal discussion of 
the major computer components. Hiring the past quarter, the following 10 groups 
totalling 49-1 people visited the computer Installation. Those attending the MIT 
summer courses represent businesses, colleges, and government agencies all over 


the United States and several other countries 

June 26 MIT Summer Course 6.569, "Switching Circuits" 18 

June 27 MIT Summer Course 2.219, "Vibration. Shock, and 94 

Noise" 

June 28 MIT Summer Course 6 569, "Switching Circuits" 21 

July 3 MIT Summer Course 16.259, "Aerodynamic Heating of 45 

Aircraft Structures In High Speed Flight" 

July 7 MIT Summer Course 15.579, "Operations Research" 33 

July 9 MIT Science Teachers Summer Program 67 

August 9 MIT Summer Course 15,548, "Electronic Computers 40 

and Business Problems" 

August 10 MIT Summer Course 1.559, "Structural Design for 52 

Dynamic Loads" 

August 15 MIT Summer Course 15.549, "Business Management 58 

and Electronic Data Processing" 

August 29 MIT Training Course for the IBM 704 66 


The procedure of holding Open House at the Laboratory on the first, and 
sometimes the third, Tuesdav of each month has continued. A total of 125 people 
attended the 4 Open House tours during the quarter, representing members and friends 
of the MIT community, Harvard, University of Virginia, Yale, New York University, 
Johns Hopkins University, Sperry Hand, Lybrand Ross Brothers and Montgomery, Shell 
Development Co., Lahev Clinic, National Shawmut Bank, Provident Institution for 
Savings, Latrola Steel Co., Ethyl Corp., Army Signal Corps, Pittsburgh Plate Glass 
Co., Dow Chemical Co., and Hawaiian Airlines. 

During the past quarter, there were also 90 individuals who made brief 
tours of the computer Installation at different times. Represented by thoso indivi¬ 
duals were Computing Devices of Canada, Rutgers University, University of Illinois, 
Stanford University, Barbour, Stockwell A Co., North American Aviation, Union 
Carbide and Carbon Co., Laval University, Quebec, Phllco Co., British Broadcasting 
System, Institute for Advanced Study, Federal Telecommunications Laboratory, 

Brandels University, Detroit Public School System. Stewart A Lloyd. Ltd , University 
College, London, Allied Control Co., Wasado University, Japan. Datamatlc Corporation, 
Metropolitan Life Insurance Co., U.S. Rubber Co., Monsanto Chemical Co., and the 
Radio Corporation of America. 
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3. ACADEMIC 


1. Introductio n 

There were a number of special summer programs offered this year which 
made use of the computing facilities In either the Office of Statistical Services 
or the Digital Computer Laboratory, respectively. The following list of subjects 
Involved direct use of the computers: 


Description 


Instructor Compu ter 


2.138 Fundamentals of Strain Gauge Techniques 

2.219 Vibration Shock and Noise 

2.769 Control Systems Engineering 

6.539 Analog-to-Digital Conversion 

6.569 switching Circuits 

8,002 Science Teachers Program 

15.548 Electronic Computers In Business Problems 

15.549 Business Management and Electronic Data 

Processing 

15.579 Operations Research 

16.259 Aerodynamic Heating of Aircraft Structures 
in High-Speed Flight 


W. Murray 
S. Crandall 
w. Seifert 
A. Susskind 
S. Caldwell 
G. Wadsworth 
R. Gregory 

R. Gregory 
P. M. Morse 


ww 

ww 

650 

WW 

WW h 650 
ww 

650 


650 

WW & 650 




In certain casei, auch aa Subject 2.769, 13.548, 15.549 and 15 579, Dr. Verzuh pre¬ 
sented a number of lectures to the classes Involved. Similarly, Professor Arden 
preaented a number of lectures in 15.579. In other caaea, the computers were used 
either on a demonstration or laboratory basis. 

2. Statistical Services Summer Academic Progra m 

There were a number of different academic projects conducted in the 
Office of Statistical Services during the summer of 1956. Specifically, these in¬ 
cluded: 

a. 6.68 - Special problems in Electrical Engineering, 

b. Individual staff member research on the development of interpretive 
systems, 

c. Master's theses research using the 650 Calculator, 

d. Foreign Student Summer Project (FSSP) Participation 


Subject 6.68 


Mr. C. W. Turk enrolled in 6.68 and carried out a 15-unlt program con¬ 
cerned with the following project: Evaluation of the M.I.T. Selective System 
(MITSS) for the Type 650 Magnetic Drum Calculator", During the course of the sum¬ 
mer, Turk evaluated the MITSS System and obtained quantitative Information regarding 
the time required to effect the following subroutines: 


a. floating addition, 

b. floating multiplication, 

c. floating division, 

d. fixed and floating square root, 




e. fixed and floating sin x, 

f. fixed and floating cos x, 

g. fixed and floating exponential, 

h. fixed and floating logarithm, 

i. fixed and floating arctangent, 

respectively. MITSS was designed to minimize storage requirements -only 405 re 
glsters are used. 


4. 650 Routines A vailab le in Stat is tica l Services 

It should be mentioned that at the present time there are a number of 
Interpretive systems available in the Office of Statistical Services for use with 
the 650 Magnetic Drum Calculator: 

a. M1TSS-MIT Selective System (405 registers), 

b. MITILAC - floating decimal Interpretive mnemonic coding (1600 regie 

ters), 

c. NOPI - fixed interpretive system numeric coding (600 registers), 

d. FLIMSY - floating Interpretive matrix system (815 registers), 

e. BTL - fixed and floating Interpretive numeric coding (1000 registers), 

f. - symbolic interpretive routine (768 registers). 

SOAP is actually a two-pass operation and does not require Internal storage during 
machine operation. In addition to the above Interpretive systems, there are a 
number of special-purpose programs which have been developed for uso on the 650 
Calculator. No attempt will be made to describe these here since a separate sta¬ 
tistical Services Report has been written to describe them. For example, REPORT 
S-29, entitled M IT Se lectiv e System for the 650 Magnetic Drum Calc ulato r. 


5. MIT Theses Using the 650 Calculator 

ftirlng the summer there were nine theses completed using the 650 Calcu 
lator and facilities of Statistical'Services. These are: 


Student 

Title 

Degree 

e upervipor 

R. E. Lyons 

Automatic Coding of a Digital 

Computer by the Compiler Technique 

Engineering 
(E.E.) 

Verzuh 

S, Matsa 

Evaluation of the Methods for Solving 
Simultaneous Linear Equations on a 
Digital Computer. 

M.S. 

(E.E.l 

Verzuh 

A. Sage, Jr, 

An Evaluation of Digital Computer 
Methods for Solving Ordinary 
Differential Equations Associated 
with a Spheroidal Antenna 

M.S. 

(E.V.) 

Verzuh 

R. Hamlin 

An Electronic Data Processing System 
for Savings Banks 

M.S. 

(Dus.Adm.) 

Verzuh 

C. H. Roth 

Study of Series-Parallel Switching 
Functions 

M.S. 

(K.E.) 

Huffman 
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C. Woodruff 

Kvaluatlon of a Method of Time 

Domain Synthesis 

M.S. 

(E.E.) 

Gulllemln 

C. Sherbrooke 

A Digital Computer Analysis for 
Production Scheduling 

B.S. 

(Bus.Adm.) 

Arden 

R. Pollett 

Solution of a Rectangular Game 

Using the 650 Calculator 

B.S. 

(Math) 

Wadeworth 

D. Dinner 

A Digital Computer Program for the 
Computation of Commutation Tfebles 

B.S. 

(Math) 

Franklin 

The actual number of machine hours required to complete a particular thesis 
varies considerably depending upon the type of thesis. Generally speaking, a 
master's thesis requires much more time than a bachelor's thesis; on the other hand 
certain theses require more machine usage because of their nature. 


6 . Foreign Student S umer ProJect 

During the summer, the following Individuals spent the summer months In 
the Office of Statistical Services: 

Dr. Heinz Cumin (SC.D.) has received his doctorate from the University of 
Muenster, Germany In 1954. He has spent the last two years working on the logical 
design for a transistorized computer for Sienens & Halske, Munich, Germany. 


August 20-31, 1956. This course was an 
of instruction over a two-week period, 
four different New England Universities 
basis This program was one of the firs 
Computation Center. 


intensive program consisting of sixty hours 
One hundred and twenty people from twenty- 
attended this program on a special invitation 
t educational contributions of the new M.l.T 


8 . M,I,T. Computation Center 

The M.l.T. Computation Center represents a cooperative venture between 
M.I.T., at least 23 New England Colleges and Universities and the I.B.M. Corpora 
tlon• The following colleges and unlverslties--ln addition to M.l.T.—are currently 
participating in the program: 


Boston College 

Bowdoin College 

Brandeis University 

Brown University 

Connecticut, University of 

ftirtmouth College 

Harvard University 

Mount Holyoke College 

New Hampshire, University of 


Northeastern University 
Khode Island, University of 
Tufts University 
Vermont, University of 
Wellesley College 
Wesleyan University 
Williams College 
Worcester Polytechnic Institute 
Yale University 


It Is expected that additional New England Colleges and Universities will become 
active participants in the very near future. 


During the summer he obtained first-hand knowledge of the applicatlonal 
aspects of computer design, using the 650 Magnetic Drum Calculator on several 
scientific problems. He Is now revising the specifications for his transistorized 
computer In light of the applicatlonal experience he received--emphasizing the im¬ 
portance of Index counters, cycling, etc. 

Mr. Ounnar Dundblad, of the Dator Corporation, Sweden, spent a number of 
weeks In the Office of Statistical Services becoming falliar with the available 
equipment. His particular Interest was that of a developmental nature, rather than 
applicatlonal. Hence, his research was devoted to transistors, printed circuit 
components, reliability studies, etc. 


The principal objective of this program Is to increase the number of 
students, engineers and scientists Qualified to use computers in their research 
efforts. The Type 704 Computer will be used solely for non-classlfled, publishable 
work originating at M.l.T., the participating New England Colleges and Universities 
and the I.B.M. Corporation, respectively. 

The 704 Computer and associated electronic accounting machine (EAM) equip 
ment will be located in a two-story annex which has been added to the new Karl T. 
Compton Laboratory (Building 26) which Is currently being erected. According to 
present plans, the building should be completed In February, 1957 and the EDPM 
equipment will be in operation shortly thereafter. The actual machine complement 
which will be available for use In the M.l.T. Computation Center Is as follows: 


Lt. DJalma S. Ferreira spent the summer In the Office of Statistical 
Services studying Evaluation of the Soots of Nonsymmetrlc M atrices . The entire 
work in this area was performed using the 650 Calculator. 4 report describing this 
work is available. 

It Is apparent from the above activities that the summer program In 
Statistical Services was extremely active and very beneficial both from the 
participants and the operational standpoint. Specifically, these studies revealed 
that a large amount of concentrated time and effort of qualified personnel Is 
needed to achieve effective solutions to these problems. A number of reports have 
been prepared describing the work of the above individuals. 


7. Specia l 704 Programming Course 

In addition to the previously-described academic programs, the first MIT 
course of programming instruction on the 704 Computer was offered during the period 


Machine Complement in the M.l.T. Computation Center 


Quantity 

Type 

Description 

i 

704 

Analytical Control Unit 

i 

711 

Punched Card Reader 

i 

716 

Alphabetic Printer 

i 

721 

Punched Card Recorder 

i 

733 

Magnetic Drum Unit (8192 words) 

i 

736 

Power Frame #1 

2 

737 

Magnetic Core Storage Units (8192 words) 

1 

740 

CRT Output Recorder 

1 

741 

Power Frame #2 

1 

746 

Power Distribution Unit 

1 

753 

Magnetic Tape Control Unit 

10 

727 

Magnetic Tape Units 

1 

780 

CRT Display Unit 
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Off-L ine Equipment 


1 

714 

Card Header 

1 

717 

Alphabetic Printer 

1 

722 

Card Punch 

2 

727 

Magnetic Tape Units 

1 

757 

Printer Control Unit 

1 

758 

Punch Control Unit 

1 

759 

Card Reader Control Unit 

1 

024 

Auxiliary Machines 

Key Punch 

5 

026 

Key Punches 

1 

056 

Verifier 

1 

066 

Printing Card Unit ) 

1 

068 

Telephone Signal Unit ) 

1 

077 

Collator 

1 

082 

Sorter 

1 

407 

Accounting Machine 

1 

519 

Reproducer 

1 

552 

Interpreter 

1 

650 

Magnetic Drum Calculator 


Liu ta 


Transceiver 


The above complement of machines represents the initial installation. As time goes 
or. and the number and type of applications change, the complement of machines will 
vary accordingly. Specifically, the two Type 737 Magnetic CORE Storage Units (8192 
words) will be replaced by a single Type 738 CORE storage (32,000 wrds) some time 
in 1958. In a similar manner, the amount of transceiver equipment, etc., will vary 
as need arises. 


PERSONNEL ON THE PROJECT 


COMMITTEE ON MACHINE METHODS OF COMPUTATION AND NUMERICAL ANALYSIS 

Faculty S uperv isors: 

Philip M Morse, Chairman 
Samuel H. Caldwell 
Herman Feshbach 
Jay W. Forrester 
James B. Keswick 
Chla Chiao Lin 

PROJECT WHIRLWIND 

Staff Members of the Scientific and Engineering Computations Group at the 
Digital Computer Laboratory: 


Physics 

Electrical Engineering 
Physics 

Electrical Engineering 
Mechanical Engineering 
Mathematics 


9. Preparation of Progra ms for the 704 Com puter 

Considerable thought has been given to the choice of methods for operating 
the 704 Computer in the M.I.T. Computation Center. As a result of a six-month study 
it has been decided that tho basic language to be used for programming will be the 
RHA»E language. SHARE, an organization of 704 users, was formally organized in 
August, 1955. Since that time, SHARE has done an outstanding Job of coordinating 
the efforts, developing new methods and distributing useful programs and materials 
to the SHARE membership. 

A more general algebraic compiler, sophisticated post-mortem routines and 
and elaborate Input-output selection routines are being written for New England 
users of the M.I.T. Computation Center. 


Frank M. Verzuh, Head 
Dean N. Arden 
Sheldon F. Best 
Marion Callaghan 
Frank C Helvlg 
Bayard Rankin 
Jack Koseman 
Arnold Siegel 
Jack Uretsky 
Murray Watkins 
Munroe R. Weinstein 
Alexander Zabludowaky 


F. M. Verzuh 
Assistant nirector 
M.I.T. Computation Center 
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